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(Nj ; Abstract 

> ; Through defining irreducible loop integrals (ILIs), a set of convene,- condition* for tfie regularised (onadrat- 

ically and logarithmically) divergent ILIs are obtained to maintain the generalized Ward identities of gauge 
£f) • invariance in non-Abelian gauge theories. The ILIs of arbitrary loop graphs can be evaluated from the cor- 
' responding Feynman loop integrals by adopting an ultraviolet (UV) divergence preserving parameter method. 



> 



X 



Overlapping UV divergences are explicitly shown to be factorizablc in the ILIs and be harmless via suitable 



C^j , subtractions. A new regularization and renormalization method is presented in the initial space-time dimension 
0^ ' of the theory. The procedure respects unitarity and causality. Of interest, the method leads to an infinity free 

o ■ 

' renormalization and meanwhile maintains the symmetry principles of the original theory except the intrinsic 
mass scale caused conformal scaling symmetry breaking and the anomaly induced symmetry breaking. Tadpole 
t-j^ , graphs of Yang-Mills gauge fields are found to play an essential role for maintaining manifest gauge invariance 
via cancellations of quadratically divergent ILIs. Quantum field theories (QFTs) regularized through the new 
method are well defined and governed by a physically meaningful characteristic energy scale (CES) M c and a 
physically interesting sliding energy scale (SES) fi s which can run from fi s ~ M c to a dynamically generated 
mass gap /j s = /j c or to (i s = in the absence of mass gap and infrared (IR) problem. For M c — > oo, the initial 
UV divergent properties of QFTs arc recovered and well-defined. In fact, the CES M c and SES at fi s = /i c 
play the role of UV and IR cutoff energy scales respectively. It is strongly indicated that the conformal scaling 
symmetry and its breaking mechanism play an important role for understanding the mass gap and quark con- 
finement. The new method is developed to be applicable for both underlying rcnormalizablc QFTs and effective 
QFTs. It also leads to a set of conjectures on mathematically interesting numbers and functional limits which 
may provide deep insights in mathematics. 

PACS numbers: PACS numbers: ll.10.-z, 11.15.-q, 11.15.Bt 
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INTRODUCTION 



All known basic forces of nature, 



Symmetry has plaid an important role in particle physics [ 
i.e., gravitational, electromagnetic, weak and strong forces, are governed by the symmetry principles. 
Three of them (electromagnetic, weak and strong forces) have turned out to be characterized by the 
gauge symmetry U(l)y x SU(2)i x SU(3) C . Thus the three fundamental gauge interactions among 
the building blocks or elementary particles (quarks and leptons) are mediated via the Abelian and 
non-Abelian Yang-Mills gauge fields^]. The real world of particles has been found to be successfully 
described by quantum field theories (QFTs) 0, 4]. In the meantime, the gauge invariance has been 
viewed as a basic principle in addition to the well-known basic principles of Lorentz invariance and 
translational invariance. QFTs have also been applied to deal with effective theories for composite 



fields and particles at low energies and a 
mechanics and condensed matter physics 



so critical phenomena (or phase transitions) in statistical 
i|. Nevertheless, QFTs cannot be defined by a straight- 
forward perturbative expansion due to the presence of ultraviolet (UV) divergences. Namely, the 
definition of Feynman diagrams in the perturbative expansion may be meaningless because of lack 
of convergence. To avoid such difficulties, one may modify the behavior of field theory at very large 
momentum, 1 or introduce so called regulators (more general speaking, the regularization quanta), so 
that all Feynman diagrams become well-defined finite quantities. Such a procedure is usually called 
regularization. The most important features required for the regularization are that the regulariza- 
tion should maintain the basic symmetry principles of the theory, such as gauge invariance, Lorentz 
invariance and translational invariance (or Poincare invariance) , and also preserve the initial but well- 
defined divergent properties of the theory. Many regularization methods have been introduced in the 



literature, such as: cut-off regularization 
regularization P], BPHZ regularization 
differential regularization 



7|], Pauli-Villars regularization 



J], 



Schwinger's proper time 
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P], BPHZ regularization |lj a- reg — 41,. 

131 ] . All the regularizations have their advantages and shortcomings. 
The cut-off regularization is the simplest one by naively setting an upper bound to the integrating 
loop momentum. This method is often used to treat QFTs in statistical mechanics or in certain low 
energy dynamical systems, where the divergent behavior of the theory plays the crucial role, and 
the Lorentz or Poincare invariance beomes unimportant and the gauge symmetry is not involved at 
all. In contrast, the method fails in applying to the QFTs for elementary particles in which Lorentz 
or Poincare invariance and Yang-Mills gauge invariance play an imporatnt role. This is because the 



1 If there exist infrared divergences, the behavior of field theory at very small momentum may also need to be modified 



2 



method destroys the principle of Lorentz and translational ( or Poincare) invariance and also the prin- 
ciple of gauge invariance for gauge theories. The spirit of Pauli-Villars's regularization is to modify the 
propagator. Its prescription is simple: replacing any propagator by a sum of propagators with large 
masses, and choosing appropriate coefficients so that the large momentum behavior becomes well con- 
trolled. More general modifications on the propagators may be introduced in Schwinger's proper time 
regularization. Whereas in Pauli-Villars's regularization, a set of regulator fields are usually introduced 
to modify the action of original theory. Though Pauli-Villars's inspired regularizations may preserve a 
large number of symmetries of the theory, there are still a class of Feynman diagrams which cannot be 
regularized by this approach. Such field theories include the non-Abelian Yang-Mills gauge theories 3] 
and also the non-linear a- model. In ref.[l|3], a higher covariant derivative Pauli-Villars regularization 
was proposed for maintaining the gauge symmetry in non-Abelian gauge theories. Nevertheless, it has 
been shown that such a regularization violates unitarity 13] and also leads to an inconsistent quantum 



chromodynamics (QCD) la] . BPHZ prescription relies on the theorem proved by Weinberg 191 ]. which 
states that the requirement for the actual convergence of the amplitude corresponding to any graph 
is that the power-counting should give a negative superficial degree of divergence for the complete 
multiple integral for the whole amplitude and also for any subintegration defined by holding any one 
or more linear combinations of the loop momenta fixed. The prescription does eliminate all superficial 
divergences and render the renormalized Feynman integrals to be convergent. The forest formula 
provides a powerful tool to construct proofs of renormalizability to all orders. In fact, the method 
is a regularization- independent subtraction scheme up to the finite part. To proceed the practical 
calculation so as to extract the finite part by means of the forest formula, one still has to adopt a 
regularization scheme. As such a subtraction process is based on expanding around an external mo- 
mentum, which modifies the structure of the Feynman integral amplitudes, thus the gauge invariance 
is potentially destroyed when applying the BPHZ subtraction scheme to non-Abelian gauge theories. 
In addition, the unitarity, locality and causality may also not rigorously hold in such a subtraction 
scheme. In the lattice regularization, both space and time are made to be the discrete variables. The 
method can preserve gauge symmetries of the theory, but the principle of Lorentz and translational 
(or Poincare) invariance is lost in this method. Though there is a great advantage that lattice regu- 
larization may be extended to the nonperturbative calculations by a numerical method, it may also, 
at the same time, lead to a disadvantage due to a very complicated perturbative calculation. The 
principle of dimensional regularization is to define all Feynman diagrams by analytic continuation in 
the space-time dimension parameter, when the space-time dimension is required to be the initial one, 
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the Feynman diagrams will recover the initial logarithmically UV divergences but the quadratically 
UV divergences are suppressed. This method provides the simplest symmetry principle preserving 
regularization in the perturbative calculations. Whereas it seems to be meaningless beyond perturba- 
tive expansion as it involves continuation of Feynman diagrams in space-time dimension parameter to 
arbitrary values. The method also fails if quantities are specific to the initial space-time dimension, 
such as 75 in four dimensions and the complete antisymmetric tensor e^ upa - as well as the case that 
scaling behavior becomes important. The cases include the chiral, topological and supersymmetric 
theories. For that a so-called dimensional reduction regularization was suggested J^J as a variant of 
dimensional regularization, in which the continuation from dimension d=4 to d=n is made by com- 
pactification. Thus the number of field components remains unchanged and the momentum integrals 
are n-dimensional. Such a prescription may cause ambiguities in the finite parts of the amplitudes 
and also in the divergent parts of high order corrections. The scheme has been applied particularly in 
supersymmetric models. In general, it seems to hold only at one loop level and is actually inconsistent 
with analytical continuation concerning 75 2l|. On the other hand, the dimensional regularization may 



not be applicable to deal with the dynamics of effective QFTs at a physically meaningful finite energy 
scale. The well-known example is the derivation of gap equation in the gauged Nambu-Jona-Lasinio 
model [l^. It was shown that the dimensional regularization cannot lead to a correct gap equa- 
tion. This is because the dimensional regularization destroys the quadratic 'divergent' (or quadratic 
'cut-off' momentum) term in the gap equation though it maintains gauge invariance. Instead, when 
adopting the simple cut-off regularization to regulate loop integrals, one can arrive at a desired gap 
equation, but the gauge invariance is destroyed by a quadratically 'divergent' (or quadratic 'cut-off' 
momentum) term and also by ambiguities associated with arbitrary routing of loop momentum due to 
lack of translational invariance. Another example is for the evaluation of the hadronic matrix elements 
of four quark operators in the kaon decays K — > tttt. If applying the dimensional regularization to the 
long-distance operator evaluation in the chiral perturbation theory, one arrives at a wrong sign for the 
leading order contributions when matching to the corresponding short-distance operator evaluation 
in perturbative QCD 22]. The reason is the same as the previous one, the dimensional regularization 
destroys the leading quadratic cut-off momentum terms which have the same sign as the leading order 
terms in perturbative QCD. Actually, the quadratic terms in the chiral perturbation theory were found 
to play a crucial role for understanding the AI = 1/2 rule and for providing a consistent prediction 
on the direct CP-violating parameter e'/e in kaon decays [221] . From these two realistic and interesting 
examples, it is not difficult to further understand the shortcoming of the dimensional regularization. 
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In general, the dimensional regularization suppresses the scaling behavior somehow via a cancellation 
between UV and IR divergences. Actually, in the dimensional regularization there is no 'divergent' 
contribution that can appear in a power of cut-off energy scale. The differential regularization was 
proposed as a method working directly on Feynman graphs in coordinate spacefill]. Where one substi- 
tutes singular expressions by derivatives of well-behaved distributions. However, in gauge theories, the 
unsatisfactory feature arises due to the arbitrariness of the different renormalization scales in the Ward 
identities among renormalized Green functions. Thus one has to adjust the scales in an appropriate 
way so as to preserve the gauge invariance. Lately, a so-called constrained differential regularization 
was introduced 



23| to fix the ambiguities of the local counterterms, so that the gauge invariance can 
be maintained in a convenient way. While from the calculational viewpoint, one loses the advantage 
of working directly in the momentum space. 

From the above points of view, one may arrive at the conclusion that up to now there exists no 
single satisfactory regularization that can be applied to all purposes in QFTs. It is then natural to ask 
whether one is able to find a regularization which can combine the advantages appearing in the above 
mentioned regularizations. To realize this purpose, the new regularization should match at least four 
criteria: 

(i) the regularization is rigorous that it can maintain the basic symmetry principles in the original 
theory, such as: gauge invariance, Lorentz invariance and translational invariance, except the anomaly 
induced symmetry breaking and the intrinsic mass scale caused conformal scaling symmetry breaking. 

(ii) the regularization is general that it can be applied to both underlying renormalizable QFTs 
(such as QCD) and effective QFTs (like the gauged Nambu-Jona-Lasinio model and chiral perturbation 
theory) . 

(iii) the regularization is also essential in the sense that it can lead to the well-defined Feynman 
diagrams with maintaining the initial but well-defined divergent properties of the theory, so that the 
regularized theory only needs to make an infinity free (finite) or infinity-controlled renormalization. 

(iv) the regularization must be simple that it can provide the practical calculations. 



The first and forth criteria are clearly understood. The second and third ones are in principle 
related, for which we may recall the following point of view: it has generally been assumed that 
the fundamental laws of nature are governed by a quantum theory of fields since the time when 
the electromagnetic and also the weak and strong interactions were found to be well described by 
the QFTs. However, this idea has been challenged from the studies of quantum gravity. As it is 
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known that the Einstein's theory of general relativity is not a perturbatively renormalizable QFT 
in the power-counting sense. Also because of divergent problem in the perturbative expansion of 
QFTs, one has raised the issue that underlying theory might not be a quantum theory of fields, it 



could be something else, for example, string or superstring[2j]. If this could turn out to be the case, 
even the quantum electrodynamics (QED) and QCD as well the quantum flavor dynamics (QFD) 
of electroweak interactions might be regarded as the effective QFTs. In general, effective QFTs are 
thought to be resulted from low energy approximations of a more fundamental theory. This may 



The theorem states that any 



become more clear from the so-called folk theorem by Weinberg 
quantum theory that at sufficiently low energy and large distances looks Lorentz invariant and satisfies 
the cluster decomposition principle will also at sufficiently low energy look like a quantum field theory. 
According to such a folk theorem, there likely exists in any case a characteristic energy scale (CES) M c 
which can be either a fundamental-like energy scale (such as the Planck scale Mp and/or the string 
scale M s in string theory) or a dynamically generated energy scale (for instance, the chiral symmetry 
breaking scale A x in chiral perturbation theory and the critical temperature for superconductivity), 
so that any effective QFTs become meaningful only at a sufficiently low energy scale in comparison 
with the CES M c . 

In general, there is no restriction on the CES M c . A particularly interesting case is to run the CES 
M c going to infinity. If in this case, the QFTs can remain to be well-defined by performing a renormal- 
ization, namely one can formally recover the initial divergent properties of QFTs in the perturbative 
expansion at M c — ► oo and consistently eliminate all infinities via an appropriate renormalization of 
coupling constants and fields in the original theory, we may then mention such kinds of QFTs as 
the underlying renormalizable QFTs. It is not difficult to show that QFTs in which there are only 
Yang-Mills gauge fields interacting with Dirac spinor fields can be viewed as underlying renormalizable 
QFTs. While the Einstein's theory of general relativity may be regarded as an effective theory at the 
sufficiently low energy scale in comparison with the CES which is in the order of magnitude of the 
Planck mass, i.e., M c ~ Mp. 

On the other hand, basing on the idea of renormalization group developed by either Wilson 
or Gell-Mann and Low 23], one should be able to deal with physical phenomena at any interesting 
energy scale by integrating out the physics at higher energy scales. This implies that one can define 
the renormalized theory at any interesting renormalization scale or the so-called sliding energy scale 
(SES) n s which is not related to masses of particles or fields and can be chosen to be at any scale of 
interest, so that the physical effects above the SES fj, s are integrated in the renormalized couplings 
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and fields. 

With these considerations, it becomes reasonable to conjecture that there should exist an alternative 
new regularization scheme that can realize the above mentioned attractive properties. Of particular, 
the regularization for which we are looking must be governed by a physically meaningful CES M c and 
also a physically interesting SES fx s , so that the laws of nature can well be described by a quantum 
theory of fields when the energy scale in the considered phenomena becomes sufficiently lower than the 
CES M c , and also the laws of nature at an interesting energy scale can well be dealt with by renormal- 
izing at the SES fi s which can be chosen to be the order of magnitude of the energy concerned in the 
considered process. It is this motivation that comes to the purpose of the present paper which may be 
organized as follows: in section II, by conceptually defining a set of so-called irreducible loop integrals 
(ILIs), we perform an explicit evaluation for the gauge field vacuum polarization graphs at one-loop 
order in the non-Abelian gauge theory with a general gauge. As a consequence, we arrive at two 
necessary conditions for the regularized quadratically and logarithmically divergent ILIs in order to 
maintain the gauge invariance of non-Abelian gauge theories (or to satisfy the so-called generalized 
Ward identities). In analogous to the quantum electrodynamics, the gauge fixing term (d^Afy 2 is 
found to be unaffected by the loop graphs. While it is unlike the usual case (such as the case in the 
dimensional regularization), the tadpole graphs of gauge fields are found in the most general case to 
play an important role for maintaining manifest gauge invariance. At the end of section, we present 
a set of consistency conditions which ensures the gauge invariance for all one loop graphs. We then 
explore in section III a new regularization method which satisfies the consistency conditions presented 
in section II for the regularized divergent ILIs. To be more clear, we begin with an explicit check on 
the two consistency conditions in the well-known cut-off regularization and dimensional regularization 
methods. It is then easily seen why the cut-off regularization destroys the gauge invariance and the 
dimensional regularization does maintain the gauge invariance. We then present an alternative new 
regularization method to achieve the desired four criteria mentioned above. It is shown in section IV 
that the regulators (or regularization quanta) are indeed decouple from the regularized ILIs. Con- 
sequently, we arrive at, through a numerical check, three conjectures on mathematically interesting 
numbers and functional limits. We then arrive at a consistent regularization scheme. In section V, 
we study in details the overlapping Feynmay integrals involved in the general two loop graphs and 
present an explicit evaluation for the corresponding ILIs by adopting the usual Feynman parameter 
method and a newly formulated UV-divergence preserving parameter method. The latter method 
ensures that all the Feynman parameter integrations contain no UV divergences. The explicit forms 
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and remarkable features of the ILIs enable us to deduce a set of key theorems which contain the fac- 
torization theorem and subtraction theorem for overlapping divergences, and the reduction theorem 
for overlapping tensor type integrals as well as the relation theorem for the tensor and scalar type 
ILIs. In section VI, we present a general prescription for the new regularization and renormalization 
method. In particular, an explicit demonstration is provided for two loop diagrams. The overlapping 
divergences are shown to be harmless after an appropriate subtraction. The subtraction process can 
be made to be analogous to the one in the dimensional regularization. All the Feynman parameter 
integrations are convergent. The scheme is seen to respects unitarity and causality. We show in sec- 
tion VII how to evaluate the ILIs of arbitrary loop graphs. After a detailed evaluation for the ILIs of 
three loop graphs, a general procedure for the evaluation of any fold ILIs is presented, which should 
be practically useful to make a realistic computation for arbitrary loop diagrams. In section VIII, 
we pay our special attention to the issues on the conformal scaling symmetry breaking, the mass gap 
genesis in Yang-Mills gauge theories and the quark confinement and deconfinement. Our conclusions 
are presented in the last section. Some useful formulae and technical details concerned in the text are 
presented in appendices. 



II. CONSISTENCY CONDITIONS FOR PRESERVING GAUGE SYMMETRY 



We start from the Lagrangian of gauge theory with Dirac spinor fields ip n (n = 1, • • • , Nf) inter- 
acting with Yang- Mill gauge field A® (a = 1, • • ■ , do) 



1 



C = MH^D^ - m)^ n - -F« u Fr (2.1) 
with 

F* v = d^Al - d v Al - ghbcA^Al (2.2) 

D^n = + igT a A a ^ n (2.3) 

Here T a are the generators of gauge group and f a b c the structure function of the gauge group with 
[T a , T b ] = ifabcT c - To quantize the gauge theory, it is necessary to fix the gauge by adding the gauge 



fixing term and introducing the corresponding Faddeev-Popov ghost term 



28J with the ghost fields r\ a . 



In the covariant gauge, the additional Lagrangian £ has been found to take the following form 

C! = -^Alf + d^ a {dX + gf abcV b AD (2.4) 
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where £ is an arbitrary parameter. Thus the modified Lagrangian is given by 



t = C + £ = Mi^D^ - m)^ n - -F« u Fr 



a. 



i 

T 

^Alf + rvWrf + gf abc v b A^) (2.5) 

281,3,13 for propagators and 



Based on this modified Lagrangian, one can derive the Feynman rules 



vertex interactions 

For simplicity, we begin with our considerations at one-loop level. To find out the gauge symmetry 
preserving conditions, it is of help to introduce a set of loop integrals which consist of the scalar type 
ones 

I -2a=f^ I7T 3. \~P)X2+a ' « = ~1, 0, 1, 2, ■ ■ ■ (2.6) 



and tensor type ones 



(2vr) 4 {k 2 -M 2 ) 



d k kpky 
2a = J (2^)4 (k 2 - M 2 ] 



3+a 



/d k k jj^kjyk pk(j , , 

(27r) 4 (k 2 - M 2 ) 4+a ' a = -1,0, 1,2, • • • (2.7) 

where the number (—2a) in the subscript labels the power counting dimension of energy momentum 
in the integrals. Here a = — 1 and a = are corresponding to the quadratically divergent integrals 
(h, h^u-) and the logarithmically divergent integrals (Iq, Jo^...). We will explicitly see below that all 
one-loop Feynman integrals of vacuum polarization Feynman diagrams can be expressed in terms of 
the above set of integrals by adopting the usual Feynman parameter method. In general, all Feynamn 
integrals from the one-particle irreducible (1PI) graphs can be evaluated into the above simple one- fold 
integrals. Note that the mass factor M 2 is regarded as an effective one which depends on the Feynman 
parameters and the external momenta pi, i.e., M 2 = M 2 (m 2 ,p 2 , " ')■ 

To be conceptually helpful, we may mention the above set of loop integrals as the one-fold irre- 
ducible loop integrals (ILIs) which are evaluated from one loop Feynman diagrams. In general, n-fold 
ILIs are evaluated from n-loop overlapping Feynman integrals of loop momenta ki (i = 1,2, • • • n) 
and are generally defined as the loop integrals in which there are no longer the overlapping factors 
{ki — kj + Pij) 2 (i 7^ j) which appear in the original overlapping Feynman integrals. It will be shown 
that any loop integrals can be evaluated into the corresponding ILIs by repeatedly using the Feynman 
parameter method (see appendix A) and a newly formulated UV-divergence preserving parameter 
method (see eq. (5.6)). 

We now evaluate the vacuum polarization diagrams of gauge fields to one-loop order. There are in 
general four non- vanishing one- loop diagrams (see Figures (l)-(4)). Their contributions to the vacuum 



polarization function are denoted by Tl^l ab (i = 1,2,3,4) respectively. The first three diagrams arise 
from the pure Yang-Mills gauge theory, their contributions to the vacuum polarization function may 
be labeled as Iljf) ab = I\^} ab + IT^) a6 + Ilffl ab . The last diagram is from the fermionic loop and its 
contribution to the vacuum polarization function is denoted by H.\fJ ab = H^} ab . The total contributions 
to the gauge field vacuum polarization function are then given by summing over all the four diagrams 

4 



t=i 



Gauge invariance requires that k^Uf u = Iif u k v = which should be true for any non-Abelian gauge 
group and valid for arbitrary fermion number Nj. This indicates that both parts H$ ab and II^ afe 
should satisfy the generalized Ward identities 

k fi Il$ ab = U^J ab k u = , (2.9) 
k ^ u (g)ab = u (g)ab k u = q ( 2J0 ) 

In terms of the one-fold ILIs, the vacuum polarization function H$ ab from fermionic loop is simply 
given by (for more details see appendix A) 

m} ah = -9 2 4N f C 2 5 ab [ dx [ 2I 2flv (m) - I 2 (m)g^ + 2x(l - x)(p 2 g^ - PflPu )I {m) } (2.11) 

Jo 

which shows that the gauge invariance is spoiled only by the quadratic divergent ILIs. It then implies 
that the regularization for which we are looking should lead to the following general condition 

I^Am) = \l?{m) (2.12) 



so that the gauge non-invariant terms caused by the quadratically divergent ILIs cancel each other 
and the vacuum polarization function H$ ab satisfies the generalized Ward identity and maintains the 
gauge invariance. Here the integrals with superscript R means the regularized ILIs. Obviously, under 
this condition the regularized vacuum polarization function II^" 6 from fermionic loop becomes gauge 
invariant and takes a simple form 



I' 1 

n S? = -9 2 ±N f C 2 5 ab (pV - PllPv ) J dx 2x(l - x)I*(m) (2.13) 

The gauge field vacuum polarization function Tx$u in the pure Yang- Mills gauge theory is much 
complicated as it receives contributions from three diagrams (for a more detailed evaluation see ap- 
pendix A). Summing over the contributions from the three diagrams and making some algebra, we 
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then obtain 



ab 



/ay 



= g^Ab [ dx{[2( 2i 2llv - l 29flu ) - Ar(3)(i - 2x)( l 2iiv - \hg„v ) } 

JO * 

+ [ (5/2 - 2x(l - x) - 3(1 - x)(l - 2s) ) p 2 g^ - (1 + 4x(l - x) ) P/4 p„ ] J 

+ Ar(3) [ ( (1 - x)(l + 2x + 8x 2 )/W 2 - x(2 + x - 2x 2 ) Vil p u ) 7 - ^ffg^hpo 

- ( (1 - 2x - 2(3 - 2x)x 2 )p 2 gP 1 gl - (3x - 4x 2 + 4x V (fl£p„ + <£p„) ) hp* } 
+ Ar(3) [ ( x + x 2 - 6x 4 + 4x 5 ) PlxPv - ( x + 2x 2 - 3x 3 - 4x 4 + 4x 5 )p 2 g^ v ] p 2 J_ 2 



+ -A 2 r(4) x(l - x) [ AX + " J>V (<?> + <fe) ] ^-2 PCT } (2.14) 



To simplify the above expression, it is useful to notice the following property of the integrals with 
respect to the Feynman parameter x 

dx x I[x(l - x)] = f dx (1 -x) J[x(l -x)} = 1/2 f dx I[x(l - x)] (2.15) 
o Jo Jo 

where the integrand /[x(l — x)] is the function of the combination x(l — x) and is invariant under the 

change of the integral parameter x — > 1 — x. With the aid of this property, one can derive more useful 

identities (see appendix A). 

To be more useful, we may formally express the tensor type ILIs in terms of the scalar type ones 

lopu = ja Iq g^v (2.16) 

I-2/j.v = 7«-2 1-2 9nv ( 2 - 17 ) 

Here and a 2 are quantities relying on the regularization and expected to be determined from the 
gauge invariant conditions. Adopting the above definitions and identities, the gauge field vacuum 
polarization function 11^; 6 is simplified to be 



ab 



g 2 Ci5 ab / dx { 2( 2I 2flu - I 2 g^ u ) + [1 + 4x(l - x)] I (p 2 g fiu - p^ 
Jo 



+ ^Ar(3) [ 1 + 6x(l - x)(a + 2) - 3a ]io {p 2 g ii u -p^Pv) 

+ ^ 2r ( 4 ) a -2 x(i - x) p 2 /_2 (p 2 g/iu - PuPv) 

o 

+ -AT(3) [ 2x(l - x)(a + 2) - 1 ] J PflP u (2.18) 
+ Ar(3) x(l-x)[ ( l + 4x{l-x))p ll p u -(l/2 + 6x(l-x))p 2 g llu }p 2 J_ 2 } 

It is seen that not only the quadratically divergent integrals in the first and last terms, but also the 
logarimically divergent term (the fifth term) can in general destroy the gauge invariance. To explicitly 
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obtain the gauge invariant conditions, we may further adopt the identity 

f dx [6x(l-x)-l] J =- / dx 2x(l -x)[ 1 - 4x(l - x) } p 2 7_ 2 (2.19) 
Jo Jo 

which is obtained by performing partial integration with respect to x and using the relevant integral 

identities presented in appendix A together with the following identity 

^ = -2(1 - 2x)p 2 7_ 2 (2.20) 

Thus the integral in the fifth term of eq.(2.18) can be rewritten as 

f dx [ 2x(l -x)(a + 2) - 1 ] 7 = -2 f dx x(l - x)[ 1 - 4x(l - x)\ p 2 7_ 2 
Jo Jo 

+ (a - 1) / dx 2x(l - x) I (2.21) 
J o 

Substituting the above results into eq.(2.18), the gauge field vacuum polarization function I\^} ah is 
further simplified to be 

n(s)° 6 = g 2 C x b ah (p 2 giiU - p^p u ) f 1 dx{[l+ 4x(l - x)} I 

JO 

+ -AT(3) [ ( l + 6sc(l-aO(ao + 2)-3ao)Io - 2x(l - x) ( 1 + 12x(l - x) ) p 2 7_ 2 ] 
+ ^A 2 r(4) a„ 2 x(l - x) p 2 /_ 2 } 

o 

+ 3 2 Ci5 a6 / dx { 2( 27 2MI/ - 7 29 ^ ) + Ar(3)^— PmPi/ x(1 - x) p 2 7_ 2 } (2.22) 
Jo 1 

Here the gauge non-invariant term is associated with the factor (ao — !)• It then implies that the new 

regularization for which we are looking should lead ao to be unit (i.e., ao = 1). Based on the above 

simplified form IL$ ab , it becomes clear that in order to maintain the gauge invariance for the gauge 

field vacuum polarization function TLp} ab in the pure Yang- Mills gauge theories, the new regularization 

must satisfy the following two conditions 

I^u = \l2 9^u (2.23) 
I^u = \lo 9v (2-24) 

where the superscript 7? denotes the regularized ILIs. The first condition for the quadratically divergent 
ILIs is the same as the one required from the gauge invariance of the gauge field vacuum polarization 
function H^J ab due to fermionic loops. The gauge invariance of the vacuum polarization function 
jjW" 6 f rom non-Abelian Yang-Mills gauge theories requires an additional condition concerning the 
logarithmically divergent ILIs. 
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So far, we have provided at one-loop order a general proof for two necessary conditions that must 
be satisfied for any regularization to maintain the gauge symmetry. Under these two conditions, the 
regularized whole gauge field vacuum polarization function can explicitly be expressed in terms 
of gauge invariant form 

u ab _ jr{g)ab jT{f)ab 
VL Ry,v — LL Rfj,u ' Rfiu 



g 2 C 1 5 ab (p 2 g^ - p^Pu) I dx { [1 + 4x(l - x)] I { 

Jo 



+ ^Ar(3) [ ( 9x(l -x)-l) I R - x(l - x) ( 1 + 12x(l -x))p 2 I R 2 } 
+ ^A 2 r(4) a_ 2 x(l-x)p 2 I R 2 } 

- g 2 4N f C 2 5 ab (p 2 g^-p^) [ dx 2x(l - x)I^{m) (2.25) 

J o 

Once applying again the identity eq.(2.19), the gauge field vacuum polarization function IT^^ can be 
re-expressed into a much simpler gauge invariant form 

= 9 2 5 ab (P 2 g»» - P^Pu) f dx { d [ 1 + 4x(l -x) + A/2 ] /<? 

J 

— JV/C7 2 &e(1 - x) I?(m) - 4CiA [ 1 - 3Aa_ 2 /16 ] x(l - x) p 2 I R 2 } (2.26) 

Before proceeding, we would like to address the following observations implied from the above 
general evaluations: First, only UV divergent ILIs destroy the generalized Ward identities (or gauge 
invariant conditions) for the vacuum polarization function of gauge fields, which strongly suggests 
the necessity of regularizing the divergent ILIs in order to make the divergent loop integrals to be 
meaningful. As expected, the UV convergent ILIs are not constrained from gauge invariant condi- 
tions. Second, any consistent regularization must satisfy the necessary conditions resulted from the 
generalized Ward identities (or gauge invariant conditions) fc^II^ = 0. It also becomes clear that 
under the two necessary conditions the gauge fixing term (d^A®) 2 in non-Abelian gauge theories is 
not affected by the loop graphs, which arrives at the same conclusion as the one in Abelian gauge 
theory. Third, the quadratically divergent terms may not necessarily be a harmful source for the 
gauge invariance once a rigorous regularization can be found to satisfy the necessary conditions for 
the regularized quadratically divergent ILIs, so that all the regularized quadratically divergent ILIs 
cancel each other. Last but not least, in contrast to the usual case in the dimensional regularization, 
in the most general case, the tadpole graphs of Yang-Mills gauge fields actually play an essential role 
for maintaining the gauge invariance, which can explicitly be seen from the above general evaluation 
on the vacuum polarization function. In fact, it is the tadpole graph that leads to the manifest gauge 
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invariant form of the vacuum polarization function even without carrying out any explicit integrations 
over the Feynman parameter x and the loop momentum k. 

Here we have only presented the explicit evaluation for the vacuum polarization function of gauge 
bosons. For all other two-point, three-point and four-point one loop Feynman graphs of fermions, 
gauge and ghost bosons arising from the gauge theory described by the Lagrangian eqs.(2.1-2.5), it 
can be shown that they all maintain the gauge invariance and satisfy the generalized Ward identities 
at one-loop level as long as the following conditions for the regularized ILIs hold 

I^u = \%v I2, (2-27) 

tflivpa = «3{^V} T 2 ( 2 - 28 ) 



for quadratically divergent ILIs, and 



Jfc = 4^ I2, (2-29) 

^Opupa = ^9{pu9pa} Iq (2.30) 



for logarithmically divergent ILIs, as well as 

1 



1-2^ = q^u 1% (2.31) 

I-2pvpa = ^9{pu9pa} 1-2 (2.32) 

for convergent ILIs. Where we have used the notation 

9{pu9pa} = 9pv9pa + 9pp9ua + 9pa9pu (2.33) 

We may mention the above conditions as the general regularization-independent consistency con- 
ditions for maintaining the gauge invariance of theories. It is easily shown that the dimensional 
regularization does satisfy the above consistency conditions. 

It is of interest to see that once the above consistency conditions for the regularized ILIs hold, the 
divergent structure of the theories can be well characterized by only two regularized scalar type ILIs 
Iq and 1^- In general, one may not need to evaluate those ILIs by any specific regularization scheme. 
As it has been turned out that in the underlying gauge theories the quadratic divergences arising 
from the vacuum polarization functions cancel each other, while the logarithmic divergences can be 
fully absorbed into the redefinitions of the coupling constants and the relevant fields. In this sense, 
one arrives at a regulator free scheme for the underlying theories. Mathematically, one only needs to 
prove that such a regularization exists and can result in the above consistency conditions among the 

14 



ILIs. Practically, we shall explicitly show in the following sections that there does exist, in addition 
to the dimensional regularization, an alternative regularization scheme that can lead to the above 
consistency conditions and in the meantime match to the four criteria stated in the introduction. 

Before ending this section, we would like to point out that the above consistency conditions hold 
exactly without any ambiguities, which differs from the so-called consistency relations imposed in 



the implicit regularization [3l|]. Where some local arbitrary counterterms parameterized by (finite) 
differences of divergent integrals arise, they may (or not) be determined by the physical conditions. 
The reason is that the consistency relations in the implicit regularization are imposed among the basic 
divergent integrals which are independent of the external momenta. Namely the mass factor in those 
basic divergent integrals remains the original mass of fields. This is explicitly seen from the following 
identity and truncation which have been used in the implicit regularization 



N 



(-iy(p* + 2 Pi -ky , (-i) N + 1 (pf + 2 Pi -k) N + 1 



(U2 _ m 2Y7+l 



[{k+pif-m?} ^ (k*-m 2 y+ 1 (k 2 - m 2 ) JV + 1 [(ifc + Pi ) 2 - m 



where P i are the external momenta and N is chosen so that the last term is finite under integration 
over k. It must be this truncation that causes the local arbitrary counterterms. 

III. SYMMETRY PRINCIPLE PRESERVING AND INFINITY FREE REGULARIZATION 

We now proceed to explore a possible new regularization which can preserve the basic symmetry 
principles of the theory and lead to an infinity free renormalization. Before doing so, it is helpful 
to briefly look at the well-known two simple regaularizations, i.e., cut-off regularization and dimen- 
sional regularization, through checking the consistency conditions resulted from the generalized Ward 
identities (or gauge invariant conditions). 

The regularized ILIs in the cut-off regularization are presented in appendix B. From the explicit 
forms overthere, one easily sees that 

f^/jtulcutoff = -^9ixv $ \cutoff - ~^il,v t ^ + 2^ J ^ 2 9twI ^ CUt °H ' < ~ 3 ' 1 ' ) 

R In _J_ ( M<2 _ .M 4 _ 1_\ 1 R 

A) iivVutoff - -^Qtiu h\cutoff + Q4^2 9 iw ( ^2 + ji4^ ~ 2(A 2 + M 2 ) 2 ~ 2J Z 9 ^ ' Cilto// ' ' 

where A 2 is the cut-off momentum. Obviously, the regularized ILIs by the simple cut-off regularization 
do not satisfy the consistency conditions. In fact, not only the quadratically divergent terms can 
destroy the gauge invariance, even the finite terms may spoil the gauge conditions in the simple 
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cut-off regularization. Whereas one may notice that the involved logarithmically divergent terms in 
the regularized ILIs do satisfy the consistency conditions. This is actually the case realized by the 
dimensional regularization. From the explicit results given in appendix B for the regularized ILIs in 
the dimensional regularization, it is not difficult to yield, by using the identity T(l + z) = z T(z), the 
desired relations 

tR I _ 1 T R\ jR I _ 1 T R\ fo q\ 

1 2 fj,v\Dim. — Ijil^v 1 2 \Dim. 5 ^0 fivlDim. — ~^Jliv 1 Q \Dira. \°-°) 

It then becomes manifest why the dimensional regularization preserves gauge symmetry. 

To see the difference between the cut-off regularization and dimensional regularization, one may 
make an expansion for e — > 0. Using the properties of the T-function 

r(|) = \ - ie + o(e), r(-i + £) = -r(|)/(i - 1) (3.4) 

with the Euler number = 0.5772 • • • , one has 

I?\ Dim . = M 2 [ - - 1E + 1 - ln(4^X 2 ) + 0(e) } (3.5) 

IolDim. = 7^-2 [---YE- H^M 2 ) + O(e) } (3.6) 
167r z e 

Comparing them to the results in the cut-off regularization, the divergences given by 1/e must be log- 
arithmic, i.e., 1/e — ► In A. This may become more clear from the vanishing integral in the dimensional 
regularization 

d D k 



k 2 



(3.7) 



It is this property that suppresses the quadratic divergence in the dimensional regularization. Actually, 
the dimensional regularization destroys all the divergences in a power of cut-off momentum because 
of the vanishing integrals 



/ 



d D k (k 2 ) 1 = / = 0, 1, ••• (3.8) 



which implies that divergences appearing in any loop integrals are expressed only in terms of the 
expansion by powers of 1/e in the dimensional regularization. 

We now turn to investigate the possible new regularization. It has been seen from the previous 
section that the divergences actually arise from loop integrals. It is the UV divergent loop integrals 
which destroy the generalized Ward identities of gauge invariance. This observation naturally mo- 
tivates us to propose a new regularization scheme which is supposed to regularize the ILIs of loop 
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graphs, i.e., the whole Feynman diagrams, instead of only the propagators carrying out in the Pauli- 
Villars inspired methods. Unlike the Pauli-Villars regularization, we are not going to introduce any 
regulator fields to modify the action of original theory. Alternative to the dimensional regularization, 
we shall not change the space-time dimension of initial theories. Our motivation is only to modify the 
very high energy behavior of loop integrals and /or the very low energy behavior of loop integrals if the 
integrals concern IR divergent. In fact, we do not well understand the physics at very short distances, 
especially the physics above the Planck scale. Thus it must not be surprised that at infinitely large 
loop momentum, the Feynman diagrams and the corresponding Feynman integrals become ill-defined 
due to the lack of convergence. They are actually meaningless because of the UV divergent integrals 
which destroy the basic gauge symmetry principle in gauge theories. Moreover, it is also a hard task 
in certain QFTs to deal with the physics at a very low energy scale when the nonperturbative effects 
of strong interactions become important at that scale. In particular, the same serious problem may be 
faced when the Feynman integrals also become ill-defined at zero momentum due to the existence of 
IR divergences. All of these considerations strongly suggest the necessity of regularization for QFTs. 

Here we are going to present a new regularization without introducing any additional Lagrangian 
formulation to modify the original theory as usually done in the Pauli-Villars inspried regularization. 
What we are doing is only to require the new regularization possessing the most important and 
necessary properties, namely: the regularization shall make the meaningless divergent ILIs to be 
well-defined and physically meaningful, and also to be applicable to both underlying (renormalizable) 
QFTs and effective QFTs. Furthermore the new regularization shall preserve the initial divergent 
behavior and meanwhile maintain the important symmetry principles in the original theory. More 
specifically we are going to work out the regularization which will ensure the regularized divergent 
ILIs satisfying the consistency conditions of gauge invariance even when keeping the quadratic 'cut- 
off' terms in the loop integrals. In general, the quadratic 'cut-off' terms will be found in the new 
regularization to be harmless for the underlying (renormalizable) QFTs but very important for the 
effective QFTs. Of important, we shall finally verify that the physical quantities are independent of any 
regulators (or regularization quanta) introduced in the regularization scheme, namely the regulators 
(or regularization quanta) will eventually be decouple from the well-defined regularized theory. 

The prescription for the new regularization is simple: supposing that the ILIs of loop graphs have 
been rotated into the four dimensional Euclidean space of momentum, we then replace in the ILIs the 
loop integrating variable k 2 and the loop integrating measure f d 4 k by the corresponding regularized 
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ones [A; 2 ]; and f[d A k] 



k 2 [A; 2 ], = k 2 + Mf , ( d A k^ f [cfikl = lim V <f / d 4 /t (3.9) 

J J N ' M ?7^ J 

In here Mf (1 = 1, • • • ) may be regarded as the mass factors of regulators (or regularization quanta). 
For I = 0, we shall have Mg = to maintain the original integrals. If the original integrals are IR 
divergent, one can set Mq = fi 2 to avoid IR divergent problem. In general, the mass factors Mf 
(1 = 1, ■ ■ ■ ) of the regulators (or regularization quanta) are introduced as unknown parameters and 
their values can be infinitely large. It is seen that the initial integral is recovered by letting the 
regulator mass factors Mf go to infinity, i.e., Mf — »■ oo (I = 1, • • • ). Where cf (I = 1,2, - ■ ■ N) are 
the corresponding coefficients and will be judiciously chosen to modify the short-distance behavior of 
loop integrals and remove the divergences. Specifically, in order to suppress the divergences at high 
orders, the coefficients cf with the initial condition Cq = 1 may be chosen to satisfy the conditions 

N 

limVcf(M, 2 r = (n = 0,l,..-), ^ = 1 (3.10) 
so that the regularized integrals of f d 4 k(k 2 ) n vanish 

N 

/ [d%i (k 2 + Mf) n = lim V cf / d 4 k (k 2 + Mf) n = (n = 0, 1, • • • ) (3.11) 
J N,Mf ~^ J 

which coincide with the ones in the dimensional regularization. Where the integer N counts the 

corresponding numbers of regulators (or regularization quanta) and is in principle a free parameter. 

Its value can be infinitely large. Practically, the regulator number iV may also only need to reach 

a sufficiently large value around which the regularized ILIs become insensitive to the choice for its 

specific values. Finally, we shall prove that the regularized ILIs must not be sensitive to the choices for 

the various mass factors Mf of regulators (or regularization quanta) and the corresponding coefficients 

cf 1 as well as the regulator number N. 

We would like to point out that the order of operations for performing the integration over the loop 

momentum and making the limitation to the regulator mass factors is very important for divergent 

integrals. One should first carry out the integration over the loop momentum, then take the infinite 



limit for the regulator mass factors, (Mf — > oo (I = 1, • • • ). Physically, this order of operations implies 
that 

M 2 

— i -> at k 2 -> oo and Mf -> oo (I = 1, 2, • • • ) (3.12) 

Namely, the loop momentum in the Euclidean space goes to infinity in a much faster way than the 
mass factors of regulators ( or regularization quanta) . 
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Before a detailed discussion, one may notice that the present regularization scheme will result in 
a quite different property for the quadratically divergent integral in comparison with the dimensional 
regularization. This is because in the dimensional regularization, the quadratically divergent integral 
vanishes for the massless case (see eq.(3.7)). While in the present regularization scheme, only when 
the loop momentum k becomes sufficiently large with k 2 3> M 2 , one then approximately has 

r°° 1 r 1 00 / M 2 \ n 

/ V2TW2 - / [*% 4 E hS- = ^ 3 - 13 ) 

which will lead to a distinguishable property in the present regularization scheme. We shall come to 
a more detailed discussion below. 

Let us first apply the new regularization scheme to the scalar type divergent ILIs 

h R = * /g| 1 . S -i hm fcf I** — L_ 
2 J (2vr) 4 -k 2 -M 2 n,m?^ J (2vr) 4 k 2 + M 2 

N 

lim E c ^ ^ ln ^ ( 3 - 14 ) 



16vr 2 JV,Mf z=0 



[<™]» 1 , ; „ V-JV f d*k 



jR = i I L" Mj 1 = f lim Vc" 



(2tt) 4 (-k 2 -M 2 ) 2 N,M?f^ Q J (2vr) 4 (k 2 + M 2 ) 2 



N 



lim Vcf InM, 2 (3.15) 

>r nj2 — 4 



167T 2 7V,M ; 2 i=0 



with 



M? = M? + M 2 



To check whether the above new method satisfies the consistency conditions for maintaining gauge 
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invariance, we further evaluate the corresponding tensor type ILIs by the same regularization scheme 

[d 4 % -k^k 

v _ . ,. n f d k k^kp 



jR = a / l_ ^ = _a lim V c N / — — M 

2 » v J {2ttY (_F _ M 2 ) 2 - NM? 1 J (2tt) 4 (Jfc2 + Mf)2 

= -0u„ o lim > / d 2 fc 1 + 1- + L 

4^" 16vr 2 N,Mf ^ Q 1 J k 2 + Mf (k 2 + M 2 ) 2 J 

1 -' ^ 

2 ^ 167T^ N,M? j—i 

1=0 

1 r* 



J, 



R 



o<fo* I? (3.16) 

TV 

lim Vcf / 



(2vr) 4 (_jfe2 _ m 2 )3 = VSf ~ 1 7 (^O 1 (fc 2 + M 2 



1 * v Ar /\ 2 , r 1 



4^" 16vr 2 iv,Mf ^ W k 2 + M 2 (k 2 + M 2 ) 2 (fc 2 + M 2 ) 3 J 

1 -i * 

-Quu n hm Vcf InM, 2 

4^ 16vr 2 JV.M'f^ ' 1 

1 1=0 

I R (3.17) 



which shows that the method does preserve gauge invariance. Note that in obtaining the above results, 
only two conditions have actually been used, i.e., 



JV 

lim <f ( M iY = with n = 0, 1 (3.18) 

1=0 



N,Mf 



To be able to directly compare with the cut-off regularization and dimensional regularization, and 
also to explicitly check whether the present new regularization maintains the harmless and important 
quadratic 'cut-off' terms, it is of help to find out the manifest expressions for the mass factors M 2 
of regulators (or regularization quanta) and the corresponding coefficients cf. For a self consistent 
reason, the explicit forms of M 2 and cf 1 are assumed to be uniquely determined through the conditions 
given in eq.(3.11). This requires that for a given N the number n of conditions should be sufficient 
in comparison with the unknown variables 2N which consist of N's mass parameters M 2 and N's 
coefficient functions c? . Also for an economic reason, it is supposed that the regularization scheme 
only involves a minimal set of parameters. With these considerations, we come to the following simple 
choice for the mass factors M 2 of regulators (or regularization quanta) 

M? = n 2 s + IM 2 R Z = 0,1,--- (3.19) 

Here Mr may be regarded as a basic mass scale for all regulators (or regularization quanta). The 
energy scale fi s may be viewed as an intrinsic mass scale and will be seen to provide an IR cut-off 
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energy scale. In the absence of IR divergences, the scale fj, s can in principle run to zero /i s = 0. In this 
case, there is only one free parameter Mr. In general, even without IR problem the scale fi s can still 
be introduced to play the role of sliding energy scale (SES). A particularly interesting case is that a 
mass gap around the scale /i s = \i c may be generated due to strong interactions, if this case happens, 
the SES [i s is naturally set to be at the mass gap fi s = fi c . 



Substituting the above explicit form of M? to the mass-coefficient conditions given in eq.(3.11) 
for the regulators (or regularization quanta), we arrive at the interesting constraints purely for the 
coefficients cf 



N 



lim V cf l n = with n = 0, 1, ■ ■ ■ N - 1 (3.20) 

7 — i.rvn ' 



1=0 



For a given N, the coefficients cf are completely determined by solving iV's linear equations with the 
initial condition cf = 1. The result is explicitly given by 

f-t-Vw^m (3 ' 21> 

which may be regarded as the sign-associated Combinations. More precisely speaking, (— l) l cf is the 
number of combinations of N regulators (or regularization quanta), taken / at a time. The initial 
condition cf = 1 is required so that we shall recover the original ILIs when no any regulator is 
introduced. 

Up to now, we have completely presented an explicit new regularization scheme in which there 

exists no free parameters when both the basic mass scale Mr and the number iV of regulators (or 

regularization quanta) are taken to be infinitely large, i.e., Mr, iV — > oo, and the SES fi s is set to 

be zero fi s = 0. What we shall do next step is to prove that the regularized ILIs are insensitive to 

the regularization scheme. More specifically the regularized ILIs should eventually be independent of 

the basic mass scale Mr and the number iV of regulators (or regularization quanta) when they are 

approaching to the infinitely large limits or practically when they are taken to be sufficiently large. To 

see that, let us apply the above explicit forms for the mass factors Mf of regulators (or regularization 

quanta) and the coefficients cf to the regularized ILIs. After some algebra, we obtain the following 

explicit and simple forms for the regularized ILIs 

? M 2 ii 2 

h R = - Y ^{M 2 c -A^ 1 f-l W + l + y2({ r p]} (3.22) 

i M 2 ii 2 



'-^-li^t 1 -^)] < 3 ' 24 > 
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where we have introduced the following definitions and notations 

N 



M -=^ E C ^ ( l ^ l ) M R ( 3 - 25 ) 

' R 1=1 

,,2 N ..2 00 ( \n-lr / „2 \ " 

yo (id) = - hm y <f ln(1 + Jd) = y id d ( id ) (3.26) 

c Z=l R n=l v c 7 

yi( id)_i = Hm l ^A Hl + J_ ) = y (-) ±2 (id) -1 (3.27) 

yu M 2; n,m r ^ 1 u y IM 2 ' (n + l)n\ \M 2 J v ; 

c j_j r K n= l \ ' St./ 

/, 2 m 2 /, 2 00 i'-'i™- 1 r / t, 2 \ n 

,, 2 - // 2 // 2 00 r "i™- 1 r / ;/ 2 \ " 

Z=l w " n=l 



fi 2 = fi 2 s + M 2 (3.30) 



Here 7^ and L n (n = 1, 2 • • • ) are the numbers of functional limits 

iV AT 



7 m = lim7 to (iV) = lim{ ^cfln/ + ln[ ^cf ZlnZ ]} (3.31) 

i=i z=i 

L n = ]imL n (N) = -lim^cf £ [ £ c ^ ^ln*' ]» (3.32) 

2=1 i'=i 



Note that the following notations are understood 

lim = lim , lim = lim (3.33) 

It is manifest that the energy scale ll s sets an IR 'cut-off' for the case M 2 = which could occur 
in a massless theory with imposing on mass-shell conditions, and the energy scale M c provides an UV 
'cut-off'. In fact, M c plays the role of characteristic energy scale (CES) and its magnitude is supposed 
to be known for the considered quantum theories, and jjL s plays the role of SES and it can in principle 
run from fi s = M c to li s = if there is no IR divergent problem and no mass gap generated in the 
considered theory. Of particular, the method does maintain the quadratic 'cut-off' term M 2 , which 
makes it different from the dimensional regularization. Though this feature is analogous to the naive 
cut-off regularization, while it is quite unlike the naive cut-off regularization as the new regularization 
ensures the consistency conditions for the regularized scalar and tensor type ILIs. So that the new 
method preserves the basic symmetry principles of the theory, such as gauge invariance and Lorentz 
invariance as well translational invariance, which is compatible to the dimensional regularization. It 
then becomes clear that the UV 'cut-off' scale M c and the IR 'cut-off' scale ll s = fi c in the new 
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regularization are completely distinguishable from the ones in the naive cut-off regularization. In the 
latter case, the cut-off momentum scales are introduced by imposing kinematically the upper and 
lower bound of the integrating loop momentum, which explicitly breaks the translational and Lorentz 
invariance and also destroys the gauge invariance in gauge theories. In the former case, which is in 
contrast to the latter one, the energy scales M c and [i s are entered intrinsically as two basic mass scales 
from the regulators (or regularization quanta). For convenience and also to be distinguishable from 
the naive momentum cut-off, the energy scales M c and jj, s may be mentioned henceforth as intrinsic 
mass scales. 

It is particularly noted that when taking the CES M c to be infinite, i.e., M c — > oo, we then recover 
the initial but well-defined UV divergences of the ILIs including the quadratically divergent terms. For 
non-abelian gauge theories of Yang-Mills fields interacting with Dirac spinor fields, as all the quadratic 
divergent terms cancel each other in the new regularization scheme (which has explicitly been shown 
in the previous section for the gauge field vacuum polarization function at one loop order), thus the 
resulting divergent behavior in the regularized gauge theories at M c — > oo is actually equivalent to the 
one by the dimensional regularization once one makes a simple replacement ln(M c //x) <-> 1/(4 — D) 
with taking D — > 4 and M c — > oo. Therefore, such non-Abelian gauge theories may be viewed, 
according to our previous convention, as underlying renormalizable QFTs. 

Before ending this section, we would like to emphasize that the new regularization described 
above differs in principle from the dimensional regularization. The differences arise not only from 
the quadratic 'cut-off' term, but actually from all the divergences and their origins. This is because in 
the dimensional regularization the divergences are intimately correlated to the space-time dimensions. 
Therefore it is inevitable to meet divergences in the dimensional regularization when one is forced to 
take the space-time dimension D to be four (D = 4) of the real world or to be the one in the original 
theories. Unlike in the dimensional regularization, the space-time dimension in the new regularization 
is the same as the one in the original theory, the initial divergences of ILIs are recovered and well- 
defined by taking the CES M c to be infinite, M c — > oo. In principle, M c is a free but well-controlled 
parameter and its value relies in general on the considered phenomena described by the quantum the- 
ories. In this sense, the CES M c is physically meaningful, so do the quadratic terms or other possible 
high order terms of M c . On the other hand, the finite terms in the new regularization are given by 
the polynomial of fi 2 /M 2 rather than the one of 1/ln ^ ~ (4 — D) in the dimensional regularization. 
Of particular, in the new regularization, one can always make on mass-shell renormalization even for 
a massless theory, like QCD with massless quarks. This is because one can always choose the SES fj, s 
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to be at an interesting finite energy scale. Nevertheless, the new regularization is formally compatible 
with the dimensional regularization as they satisfy the same consistency conditions. Finally, we want 
to point out that the prescription of the new regularization formally appears to be similar to the one 
of Pauli-Villars, but they are different in the following two aspects: First, the new regularization is 
applied to the ILIs of Feynman loop graphs rather than the propagators imposed in the Pauli-Villars 
regularization which manifestly spoils gauge invariance. This is a conceptual extension to the Pauli- 
Villars regularization. Second, in the new regularization, the coefficients cf (I = 1, 2, • • • ) are well 
determined and independent of the regulator mass factors Mf for Mf = p% + IM 2 R (I = 1,2,...). 
Here the basic regulator mass scale M R and the number N of regulators are taken to be infinitely 
large in the new regularization, which is basically a crucial generalization to the usual Pauli-Villars 
regularization in which the regulator mass scale M R is given Mr = A and the number of regulators is 
normally chosen to be only a few limited one, say N = 2. As we will show below that for a finite N and 
given M R , the integrals are regulator dependent (one can explicitly see from table 1 in appendix C). 
Only when taking the number A" and the mass scale M R of the regulators to be infinite (or sufficiently 
large), then the resulting regularized theory becomes independent of (or insensitive to) the regulators. 



IV. DECOUPLE OF REGULATORS AND MATHEMATICALLY INTERESTING NUM- 
BERS AND FUNCTIONAL LIMITS 

One may notice that the basic mass scale Mr of the regulators (or regularization quanta) does 
not appear in the final expressions of the regularized ILIs, while it is intrinsically correlated to the 
regulator number A" in such a way that: taking Mr and N to be infinite, Mr — > oo and A" — > oo, but 
keeping the CES M c remain held fixed. This implies that Mr must approach to infinity in terms of 
an appropriate function of N. To obtain its explicit form, we begin with the conjecture that 

M| = M 2 In AT (4.1) 

Here M a is supposed to be a constant mass scale. Substituting this form into the definition of the 
CES M c , we arrive at the relation 

N N 

M c = E c ^ ( z ln 1 ) m r = M o % 0* N ) E c ^ ^ ln o = M o/ h ™ ( 4 - 2 ) 
' r i=i i=i 



with h w being given by 



h w = lim h w (N) = lim — jr, (4.3) 

at^oo wy ' N^°°Y,i=i<f (nnZ)lnAT ^ ; 
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Here the function h w (N) is purely determined by the regulator number N and supposed to have a 
finite limit at N — > oo. 

So far the regularized ILIs are completely characterized by the three type of functional quantities 
h w (N), 'yw(N) and L n (N) (n = 1,2, •••) , they are solely determined by the regulator number 
N. It then becomes clear that in order to prove the regularized ILIs to be independent of the 
regulators (or regularization quanta), one actually only needs to show that the regularized ILIs 
are not sensitive to the regulator number N at sufficiently large N. This is then equivalent to 
check whether the functional quantities h w (N), 7«,(iV) and L n {N) approach to certain finite limits 
at N —¥ oo, namely whether h w , ^y w and L n are eventually reached to the constant numbers. To 
be manifest, we perform, instead of an analytic proof, a numerical evaluation. The numerical 
analysis is presented in appendix C . The numerical results given in table 1 of appendix C truly 
show that for a sufficiently large TV the functions h w (N), j w (N) and L n {N) do approach to 
certain values which become less sensitive to the changes of the number N once N is going to be 
sufficiently large. In general, we are led to the following three conjectures which have been checked by 
an appropriate functional fitting method up to the precision of 1(T 3 (for more details, see Appendix C). 



Conjecture 1. 



h w = lim h w (N) = 1 (4.4) 

N— too 



which implies the following interesting functional limit 

N 



W N ee J2(-l) l {N m iy u (Ihxl) N =°° 1/lniV (4.5) 



Conjecture 2. 

lw = lim lw (N)= lim (B-^TF^nwy lnl-Hh w (N)lnN 



Nl 



= 0.5772- •• = lE (4.6) 

Namely ^y w is equal to the Euler number 7^ = 0.577215664901533 • • • . It then implies another inter- 
esting functional limit 

N m 

{N _' iyn ln/ N =°°lnlniV + 7i , (4.7) 



Conjecture 3. 



L n = lim L n (N) = 1 , n = 1, 2 • • ■ (4.8) 

N— >oo 
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which indicate the functional limits 

p(») _ Vr iV-i JV! 1 N -°° ( lniV )" „ - 1 2 C4 

^ =Z.(- 1 ) (AT _/),/. P " n-1,2- (4.9) 

When the above conjectures hold, the functions yi{x) defined in the previous section via the poly- 
nomials can now be expressed by the following simple functions 

y_ 2 (x) = l-e~ x (4.10) 

rx 1 _ -a 

y (x) = / da (4.11) 



o 



yi(x) = - {e~ x - 1 + x) (4.12) 



1 

x 

j/2 (z) = yo(^) - yi(x) (4.13) 



It is then clear that the regulators (or regularization quanta) are indeed decouple from the regular- 
ized ILIs. Furthermore, the consistency conditions eqs.(2.27-2.33) do hold in the new regularization 
(one can easily check that by using the useful integrals given in appendix D). Therefore it can be 
concluded that we do arrive at a consistent description for the new regularization scheme. 

Obviously, an analytical proof for the above conjectures must be very helpful and important. It 
may also provide deep insights in mathematics. 



V. FACTORIZATION OF OVERLAPPING DIVERGENCES AND 
REDUCTION OF OVERLAPPING TENSOR INTEGRALS 

The gauge symmetry preserving and infinity free regularization described above is well applied to 
regularize one loop graphs and the relevant ILIs. In order to obtain a consistently regularized theory, 
such a regularization should be applicable to all loop graphs. Furthermore, though the regularized 
loop graphs become well-defined and infinity free, the coupling constants and the fields in the quantum 
theories still need to be renormalized. The only difference is that the subtraction terms become finite 
quantities instead of infinite ones. The procedure of such subtractions (or the so-called renormaliza- 
tion) must be consistently carried out order by order in perturbative expansion. In fact, whether the 
subtraction terms in the renormalization are finite or infinite, in a given order any new subtraction 
terms must be finite polynomials in external momenta in order to maintain unitarity of the theory. 
Such a statement should in general be independent of any regularization. Therefore, one shall be able 
to prove this statement in a more general way without involving any regularization. 

To be specific and also for an explicit demonstration, we will treat in this section the Feynman 
integrals involved in two loop diagrams. The generalization to more loop graphs is straightforward. 
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It is different from one loop graphs, two loop diagrams require us to make a careful treatment for 
overlapping divergences. Eventually, we must show that after appropriate subtractions of one loop 
divergences, the coefficient functions of all the divergent terms are finite polynomials in the external 
momenta. Namely one has to prove the theorem that all eventual divergences must be harmless 
divergences. To arrive at this purpose, we shall first prove that overlapping divergences can completely 
be factorized. More specifically we shall show that the overall divergences of two loop graphs and the 
sub-integral divergences corresponding to one loop graphs can be well separated. It will be seen 
below that all of these properties will become manifest after performing the evaluation of two-fold 
ILIs involved in two loop diagrams. 

The explicit demonstration in the previous section shows that all divergences in one loop diagrams 
are harmless ones. In fact, as all one loop integrals can be expressed in terms of the one- fold ILIs 
by using the Feynman parameter method, therefore the divergences are completely characterized by 
the one-fold ILIs. As the coefficient functions of all divergent one-fold ILIs are finite polynomials of 
external momenta and masses, therefore all the divergences at one loop level are truly harmless ones. 

We now come to study two loop graphs. Consider first the scalar type overlapping Feynman 



integrals involved in the general two loop diagrams (i.e., the so-called general a/37 diagrams 111 ]) 



'Sv = / d 4 h I d% — — - — * — (5.1) 



The above integral has been given in the Euclidean space of momentum through an analytic Wick 
rotation. Here a, /3, 7 > since if one of them is zero, the integral is no longer an overlapping one 
and it actually becomes a factorizable one. In the gauge theories, the usual power counting rule shows 
that the divergent overlapping integrals of two loop diagrams are at most quadratically divergent, i.e., 
ot, P, 7 > 1. The above form of integrals can be regarded as the most general one for two loop graphs. 
Here we have supposed that one can always apply the Feynman parameter method to all momentum 
factors k\, so that all propagators concerning k\ and external momenta pf (i = 1, 2, • • • ) are combined 
into a single one. Similarly for the propagators concerning {k\ — k 2 ) and &2- Therefore, the three 
mass factors M\ , M. 2 , an d -M12 are in general the functions of masses m? and external momenta pj 
(t = l,2,..-) 

M 2 2 = M 2 2 (m 2 ,p 2 ,mlp 2 ,---), M 2 = M 2 (m 2 ,p 2 , m 2 ,p 2 , • • • ) (i = 1,2) (5.2) 
Based on the usual power counting rule, we may first recall some useful definitions following the 



ref. 



Ill ]: (i) the sub-integral over k\ is said to be convergent or divergent according to a + 7 > 2 



or a + 7 < 2, similarly, the sub-integral over k 2 is said to be convergent or divergent according to 
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/? + 7 > 2 or /? + 7 < 2; (ii) the overall integral of a(3j diagram is said to be overall convergent 
or overall divergent according to a + /3 + 7 > 4 or a + /9 + 7 < 4; (iii) a harmless divergence is a 
divergence with its coefficient functions a polynomial of finite order in the external momenta. In the 
gauge theories, the overall divergence of a nontrivial overlapping integral is at most quadratic, thus 
a + P + j > 3. 

In two loop graphs, we shall also involve the general tensor type Feynman integrals 

j(2) r d , f d4 (h,k lu , kl ,k 2 k 2 ,k 2u , ■■■) 

°^ 7 y (J^ + A<?)«(A| + ^((ibi-fc2+p) 2 + Al? 3 ) 7 

Similar to the one loop case, by adopting appropriate parameter methods, the tensor type Feynman 
integrals can eventually be expressed in terms of the scalar type integrals through the aid of the 

metric tensor g^ v and the external momentum vector p^. We shall first focus on the general scalar 

(2) 

type Feynman integrals I^L, and come back to the tensor type Feynman integrals at the end of this 
section. 

Though the above general forms of Feynman integrals may be more concise than the original 
ones due to the uses of Feynman parameter method, they are not yet the desired ILIs because of 
the overlapping momentum factor (k\ — k 2 + p) 2 - To obtain the ILIs of two loop graphs, according 
to the definition of ILIs, we shall further adopt appropriate parameter methods to treat the factor 
(ki — k 2 +p) 2 . As the first step, we still apply the usual Feynman parameter method to the sub-integral 
over k\ (one can also consider first the sub-integral over k 2 as the two sub-integrals are symmetric 
with the exchanges k\ <-» k 2 , a <-» (3, M\ <-> M\ and p <-» —p, for convention, we will henceforth make 
the treatments according to the order k± , k 2 , • ■ ■)■ 

Taking the Feynman parameter variables a\ = kf + A4f and a 2 = (k\ — k 2 + p) 2 + M\ 2 -, the 
sub-integral over k\ becomes a quadratic (see Appendix A) 

m = E£+± f l ix I Al ( a 2 *>-(!-*)- 



with definition 

Ml = M\ + x{M\ 2 - M\) (5.5) 

Where we have made a momentum shift fei^ — > ki^ + x(k 2fJi — p^). 

To further combine the denominators into a single quadratic, it is found to be of help to adopt an 
alternative identity 

1 T{a + (3) r A ^ 

du \„ , ^.i^-i-fl ( 5 - 6 ) 



i a b? r(a)rG9) J [a + bu] a +P 
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The advantage of this identity is that once u approaches to infinite in a way similar to a, i.e., u ~ 
a — > oo, the divergent structure of the integral coincides with the one of the term a. For convenience, 
we may mention this parameter method as an UV-divergence preserving parameter method, and 
the parameter u as the corresponding UV-divergence preserving integrating variable. Its important 
consequence is that by appropriately applying this UV-preserving parameter method to the combined 
propagators obtained by carefully using the usual Feynman parameter method, one can factor out all 
the UV divergences, so that the usual Feynman parameter integrations contain no UV divergences. 
Only IR divergences could be hidden in the Feynman parameter integrations. This is crucial point 
for analyzing the overlapping divergences. In fact, it will be shown next section that both UV and 
IR divergences can be simultaneously regularized in the new regularization. Thus all the Feynman 
parameter integrals are convergent. 

Taking a = x(l — x)(k<2 — p) 2 + k\ + M 2 and b = fe| + M\ aR d making a momentum shift 
&2/x — ► k2/j, + x(l — x)pfj,/{u + x(l — x)), we then obtain from eqs.(5.4) and (5.6) that 

■lr r roc 



~ T ^ + ^ + ^ [ dx f d % [ d * k9 f 



du 

x"i- 1 (l-x) a - 1 vP- 1 



[ k\ + Ml + (u + x(l - *)) k\ + u{M\ + ^p 2 ) ]«+/>-Hr 



(5.7) 



It is clear that once u approaches to infinite in a way similar to kf , i.e., u ~ k\ — > oo, the usual power 
counting rule implies that the resulting divergent structure does coincide with the one in the original 
sub-integral over k\. In other word, the UV divergent property of the sub-integral over k\ will be 
characterized by the one of the integral over u. It will become more clear below from the explicit form 
of ILIs. 

In general, a + (3 + 7 > 2, the sub-integral over k\ becomes convergent and can be carried out 
safely. It is straightforward to perform the integration over k\ , which leads to the result 

where we have factored out the coefficient factor [it + x(l — x)] of the momentum square k\ and 
introduced the definitions 

T{a + f3 + 1 ) 1 r(q + /? + 7 -2) 

a ^ T(a)T(f5)T(l) (a + /3 + 7 - !)(« + /? + 7 - 2) T(a)T((3)T( 7 ) { ' > 

»l = -—k V ( Ml - x(l - x)( Ml-p 2 )- f} 1 - X) \ p 2 ) (5.10) 



u + x(l - x) \ x u + x{l 
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It is useful to change the dimensionless integrating variable u into the dimensionful one k 2 = q 2 u. 
Here q 2 with dimension of mass square is introduced as an universal energy scale. In general, it only 
requires — q 2 < 4/i 2 in order to respect unitarity. Typically, one may take q 2 ~ fi 2 to reduce unnecessary 
parameters. The variable k 2 may be regarded as a momentum-like one in connecting with the initial 
loop momentum k\. Though kf cannot wholly be compatible with the initial momentum k 2 , whereas 
the integral over k 2 does maintain the one-loop divergent behavior of the initial loop momentum k\ . 
Using the new integrating variable and noticing the formal relation of the integrating measure 

d 4 ki 1 
/o J kf 



We can rewrite the above integral eq.(5.8) into the following form 

• 1 0-1 



Jo i=0 



d 4 h i r, 4l i 

a k 2 



k\ (k 2 + x(l - x)q 2 )«+7-l+i J \k 2 +M 2 + M | f Y+^~ 2 
i f>-l 



= T a ^ [ dxx-r-'il-xr-' (g 2 )^- 2 ^cf- 1 (x(l-x) go 2 ) 

I 



(M 2 + /4) (5.12) 



where I^^iM^ + M|a) ^ s an overall one- fold ILI involved in two loop graphs 



with 



- w^h (5 ' 14) 

4 = 15 K—Tl (mI-x(1-x)(M1-p')- f {l ; l)2 f p 2 ) (5.15) 

fc i fcf + x(l - x)q 2 \ kf + x(l-x)q 2 I 



According to the definition of ILIs, the above resulting loop integrals may be called as the scalar 
in£ 

in the infinite limit of k 2 , i.e., 



type two-fold ILIs. Basing on such two-fold ILIs and noticing the interesting property that //| 2 vanishes 

k 1 



i| 2 ^0 at k\ -> oo (5.16) 

fc i 



we arrive at the following important observations: 
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(i) For the simple case that a + 7 = 2 and (3 = 1, it becomes manifest that the UV divergent 
behavior of the two-fold ILIs at kf — > 00 and k\ — > 00 coincides with the one of the original integral. 
This explicitly shows that the one-loop UV-divergence preserving sub-integral over the momentum-like 
variable k\ does characterize the high energy behavior of the sub-integral over the loop momentum k\ 
in the initial integral. 

(ii) The sub-integral over the loop momentum ki turns out to be an overall one-fold ILI which is 
correlated to the sub-integral over the momentum-like variable only via the /c 2 -dependent mass factor 
Ml + n 2 %l . 

(iii) When kf — > 00, the two sub-integrals over £?2 and k± become factorized ones due to the property 
M?2 — > at kf — > 00. Thus the divergent properties of the initial two loop integral can completely be 
described by the ones of the divergence factorizable two-fold ILIs. 

(iv) The most divergent behavior of the sub-integral over k\ is governed by the power counting of 
a + 7 which reflects the divergence of one-loop sub-diagram, and the divergent property of the sub- 
integral over k2 in the ILIs is solely characterized by the power counting of a + j3 + 7 which actually 
describes the overall divergence of the general a{3~f two loop diagrams. 

(v) The external momentum dependence is contained only in the mass factor J\A\ + °f the sub- 
integral over ki- At the divergent point of the sub-integral over k\ (k\ — > 00), the possible external 
momentum dependence is solely given by the mass function M% which appears only in the sub-integral 
over k2, therefore such a harmful divergence is expected to be eliminated by one subtraction term. 

From the above observations, it becomes clear that in the two-fold ILIs the overall divergences 
of two loop graphs are completely factorized from the sub- integral divergences of one loop graphs. 
Consequently, we are able to deduce the following theorems without involving any regularization. 

Theorem I (Factorization Theorem for Overlapping divergences). Overlapping divergences which 
contain divergences of sub-integrals and overall divergences in the general Feynman loop integrals 
become completely factorizable in the corresponding ILIs. 

This theorem is the crucial one to treat the problem of overlapping divergences involved in two and 
more loop diagrams. With the observations (iv), (iii) and (ii) in the ILIs, its proof becomes manifest. 
It is actually a direct consequence of the ILIs. Therefore, the whole demonstration of the theorem is 
equivalent to the evaluation of ILIs from the general overlapping Feynman integrals of loop graphs. 
More explicit demonstrations will be presented in next section from the regularized ILIs. 

Theorem II (Subtraction Theorem for Overlapping divergences). The difference of the general 
Feynman loop integral with the subtraction term corresponding to the divergent sub-integral contains 
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only harmless divergences. 

This theorem is the central one to obtain a consist theory. Obviously, the factorization theorem 
(theorem I) becomes crucial to yield this theorem. Its proof is straightforward from the observations 
(i)-(v) in the corresponding ILIs of the general Feynman loop integrals. Let us present the simple 
demonstration based on the observations from the two- fold ILIs. Note that as the external momentum 
dependence in the two-fold ILIs only appears in the mass factor + /^fa °f the overall one-fold ILI 
over the loop momentum k 2 , and also as the sub-integral over the momentum-like variable k\ has 
the same one- loop UV divergent structure as the one over the loop momentum fci, one only needs to 
introduce one subtraction term of the sub-integral over the loop momentum k\ 

I {2 ^ = [ d^k m 1 Q , - ! d 4 k 2 3 , (5.17) 

where the superscript S denotes the subtraction term. With particularly noticing the property /x| 2 ~ 
0(^2 ) — > at k\ — > co, it then becomes clear that the difference of the integrals, i.e., ( — ), 
contains only harmless divergences. Its manifest demonstration will be given in next section for the 
regularized ILIs. 

Theorem III (Harmless Divergence Theorem). If the general loop integral contains no divergent 
sub-integrals, then it contains only a harmless single divergence arising from the overall divergence. 

This theorem may be deduced from the Theorems I and II. Its proof is also obvious with the 
observations (iv), (iii) and (ii) in the ILIs. Specifically, as a + 7 > 2 and (3 + 7 > 2 the integral 
over the momentum-like variable k\ is convergent in the two-fold ILIs. Only the integral over the 
loop momentum k 2 may contain divergence. As the integral over the loop momentum k 2 characterizes 
the overall divergence of the initial loop integral and it is actually an overall one- fold ILI, also as 
the sub-integral over the momentum-like variable k\ contains no external momentum-dependence, we 
then come to the statement in the theorem. 

Theorem IV (Trivial Convergence Theorem). If the general loop integral contains no overall 
divergence and also no divergent sub- integrals, then it is convergent. 

This theorem is really trivial and it is presented only for completeness. 

We now turn to the tensor type Feynman integrals. It is not difficult to arrive at the following 
theorems: 

Theorem V (Reduction Theorem for Overlapping Tensor Type Integrals). The general overlapping 
tensor type Feynman integrals of arbitrary loop graphs are eventually characterized by the overall one- 
fold tensor type ILIs of the corresponding loop graphs. 
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This theorem is the key theorem for the generalization of treatments and also for the prescriptions 
from one loop graphs to arbitrary loop graphs. From this theorem, it is not difficult to deduce the 
following theorem 

Theorem VI (Relation Theorem for tensor and scalar type ILIs). For any fold tensor and scalar 
type ILIs, as long as their power counting dimension of the integrating loop momentum are the same, 
then the relations between the tensor and scalar type ILIs are also the same and independent of the 
fold number of ILIs. 

This theorem is crucial to extend the consistency conditions of gauge invariance from divergent 
one- fold ILIs (or one loop graphs ) to more fold ILIs (or more closed loop graphs). 

In here we provide an explicit proof for the two loop case. The extension of the proof to more 
closed loops is obvious. Repeating the evaluation similar to the one for the scalar type two-fold ILIs, 
one can easily observe that all the tensor type Feynman integrals with structures such as kx^kxy • • • , 
k\n,k2v • • • and k^akiv • • ■ will be evaluated into the non-trivial tensor type ILIs with tensor structures 
given only by k 2 ^k 2u ■ ■ ■ plus the tensor type integrals which are constituted from the scalar type ILIs 
and the metric tensor g^ u as well as the external momentum vector p^. This can easily be seen from 
the fact that in evaluating the ILIs all involving operations only concern the simple loop momentum 
shifts 

fcl/i -> fci/x + x(k 2fl - Pp) k 2fl ->■ k 2fl + x(l - x)p ll / (u + x(l - x)) 

together with the integration over the loop momentum k\. Here the integration is well-defined as the 
sub-integral over k\ becomes convergent after applying for the UV-divergence preserving parameter 
method. 

We finally find that the non-trivial tensor type ILIs in two loop graphs are in general given by 

J ° i=0 

d 4 h 1 I 4 k 2fl k 2u ■ ■ ■ 

kj (kl+ x(l - x)ql )«+7-i-H J \k 2 2 +M 2 2 + )«+/3+7-2 

13-1 



= T Q/37 I 1 dxxl-^l-x)"- 1 f{x) {q 2 Y^- 2 Y,4~\^-*)<l 2 oj 
Jo <=o 

/lf (t; + ,(i-^)^- +i ^^- (>li+ ^> < 5 ' 18 > 

where f(x) is a polynomial function of Feynman parameter x and its form only relies on the considered 
tensor structure. Here the tensor type ILIs ^...(A^l + M^) are defined as a non-trivial overall 
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one- fold tensor type ILIs with an effective mass factor M\ + fi 



& ^(Mi + = y ^ r7? _M| nro (5.19) 



The superscript (1) means one-fold. It is obvious that 

j(2) 7 (1) 
[<5] /lu- - [S] \w 



r(2) r(l) 



(5.20) 



where [5] labels the power counting dimension of momentum. It can be generalized to any fold ILIs 

j(n) 

r(«) r(l) ^ ' ' 

2 [S\ [S] 

Here the superscript (n) labels the fold number of ILIs and its value is arbitrary. 

The above explicit forms of the tensor type ILIs provide the desired results and complete the proof 
for the theorem V and theorem VI. It is clearly seen that we eventually only need to consider the 
non-trivial overall one-fold tensor type ILIs. 

All of the above observations and the deduced theorems imply the importance of evaluating the 
ILIs of loop graphs. 



VI. REGULARIZATION, RENORMALIZATION AND UNITARITY 

It has been seen in the previous section that the problem of overlapping divergences is not relevant 
to any regularization. In the section II, it has been shown at one-loop level that the generalized 
Ward identities (or gauge invariant conditions) are in general spoiled only by the divergent integrals. 
This implies meaningless of divergent integrals and suggests the necessity of regularizing the divergent 
integrals. Consequently, consistency conditions between the regularized scalar and tensor type ILIs 
have been resulted from the generalized Ward identities at one-loop level. In section III, it has been 
demonstrated that the consistency conditions can provide stringent constraints on the regularization 
methods. In this section, we shall apply, as a practical computation and also an explicit check, those 
theorems obtained in the previous section to the regularized ILIs. In fact, by explicitly carrying out 
the relevant integrations in the regularized ILIs, we shall arrive at an independent verification on those 
theorems. 

We shall first generalize the regularization described in the section III for the one-fold ILIs to 
the two- and n-fold ILIs with n being arbitrary. The general prescription for the new regularization 
method is simple: 
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(i) Analytically rotate the four dimensional Minkowski space into the four dimensional Euclidean 
space of momentum by using Wick rotation. 

(ii) Appropriately evaluate Feynman integrals of loop graphs into the corresponding ILIs by adopt- 
ing the usual Feynman parameter method and the newly formulated UV-divergence preserving pa- 
rameter method. 

(hi) Universally replace in the ILIs the loop momentum square k 2 and the corresponding loop 
integrating measure f d 4 k as well as the UV-divergence preserving momentum- like variable k (or 
the UV-divergence preserving integrating parameter u ) and the corresponding integrating measure 
J d k/k (or J du) by the regularizing ones [A: 2 ]; and J[d 4 k]i as well as [k]i (or [u]i ) and f[d 4 k/k ]i 
(or f[du]i ), i.e., 

N 



k 2 -> [k 2 ]i = k 2 + Mf = k 2 + fi 2 + IM 2 R , f d 4 k^ f [d 4 % = lim ^ cf ! d A k , (6.1) 

J J N,Mr l=Q J 

m, S P + Mf = e +li + Ml , f**~ S Jta Ecf / ** (6.2) 



or 



1=0 

N 



[„], = u + M 2 1 q 2 = u+ ^ 2 + lM 2 R )/q 2 , f du - f [du] t = lim ^cf f 

J J > R 1=0 



du 



where the coefficient function q is given by eq. (3.21) and the mass scale M R is characterized 
via eqs.(4.1-4.4). The above prescription should be applicable to any fold ILIs. Here one needs 
to distinguish the difference of the integrating measures between the loop momentum k and the 
momentum- like variable k. In general, the SES fi s is taken to be a hnite energy scale in order to avoid 
possible IR problem, so that one can always make on mass-shell renormalization even for a massless 
theory. 

With the above prescription, the two- fold ILIs in eq. (5.18) are simply regularized into the following 
integrals 

r (2)R _ r I , -7-1 M _ ^a- 1 (^2^+7-2 



J$? = r ^7 f Q dx x1~\l - X)*- 1 (q 2 ) 



d 4 h,_ ^ J3-1 _(x(l-x)q 2 Y 

1 



k i 7^o 1 ( k i+ < l ~ x )<& + M i ) a+ ^~ 1+i 

[d ik 2]l' r ,o | . ,o , 7J2 r ..2 ( 6 ' 3 ) 



( k l + m 2 + Mf + »\ l+Mf )°+f*r~* 

where the superscript R means regularized ILIs. 

Based on the factorization theorem (Theorem I) and the subtraction theorem (Theorem II) for 
overlapping divergences as well as the harmless divergence theorem (theorem III) and the trivial 
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convergence theorem (theorem IV) presented in the previous section, we shall be able to decompose 
the regularized ILIs into the following general form 

T (2)R _ T (X) RD Al)RD T (1)RC A1)RD T (2)RC (R 

where the superscripts l RD J and l RC represent the regularized divergent and convergent ILIs respec- 
tively, and the numbers (1) and (2) in the superscripts label the one-fold and two-fold ILIs respectively. 
The possible harmful divergences must appear only in the first term with double divergences in the 
two- fold ILIs. 

To explicitly check the above decomposition, we may first carry out the integration over the loop 
momentum &2 as it is actually an overall one-fold ILIs. One can directly read off the results from 
Appendix D for the regularized one-fold ILIs. 

We consider first the case in which the overall integral is logarithmically divergent, i.e., a+fi+j = 4. 
The integration over k% results in the following explicit form 

Jo 

d A k, ^ b-1 {x(l-x)qly 



k l ( kf + x(l — x)q 2 + M 2 )°+7-l+* 

M 2 u 2 

{ [ In - lw + yA) ] - [ ln( 1 + ^ +M ,/^ M ) - «b(*? + M 2 ) } } , (6.5) 

We consider next the case in which the overall integral is quadratically divergent, i.e., a + /?+7 = 3. 
The integration over k2 leads to the result 



I? )R = r Q/37 vr 2 f dxx^\l-x) a - 1 (q 
Jo 

0-1 - x)q 2 ) 



J K l i=0 



0-1 \J2M 

I 

M 2 



( k\ + x(l - x)q 2 + Mf )<*+7-i+* 



Here the functions an d %i are defined as 



+ Mj ) = y ( j^-) - yo(^) (6.7) 

z Sl + M?)=y / M+ ^ M ? )-y 2{ ^) (6.8) 

Mm = A*f + A<2 (6-9) 
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From the explicit forms of the functions yo(x) and y2(x) (see eqs.(3.26) and (3.28)), it is easily seen 
that in the high energy limit the functions zq and Z2 are approaching to zero in terms of the inverse 
powers of k\ , i.e., 

M" 2 

zoikl + Mf) ~ 0{-^p-) ~ 0{j^—) - at fc? -» oo (6.10) 
g 2 (fcjf + M^)~0( i : J "' )~0( , ' f0 J ->0 at fe^oo (6.11) 



Mm ' ""kl + Mf 



It will be seen that the above properties are crucial for the treatment of overlapping divergences. 

To be more clear, consider first the simple case with a + 7 = 2 and (3 = 2. Applying the usual 
power counting rule to the general overlapping loop integrals, one sees that the overall integral in 
this case is logarithmically divergent (a + j3 + 7 = 4), the sub-integral over k\ is also logarithmically 
divergent (a + 7 = 2), and the sub- integral over ki is superficially convergent (j3 + 7 = 3). Turning 
to the corresponding ILIs, the usual power counting rule shows that the most divergent part of the 
sub-integral over the momentum-like variable k± remains logarithmically one (a + 7 = 2), whereas the 
sub-integral over k% in the corresponding ILIs becomes logarithmically divergent one (a+ f3+j — 2 = 2) 
as it characterizes the overall divergence of the overlapping loop integrals. One can easily carry out 
the integration for the regularized divergent part in the integral eq.(6.5). Indeed, we find that the 
regularized ILIs in this case can be written into the following form 

T (2)R _ f{X)RD M)RD f(l)RC M)RD J2)RC , -v 

J — J 1 1 (O.lZj 

where the involving regularized one-fold ILIs are given by 

-1 ,,2 



I$ )RD = - 2 / dx [ In ^ - lw + yo fe) ] (6.13) 

JO M 9o M c 

I^ RD = 7r2[ ln^-7„ + yo(§)] (6.14) 
Mm lvl c 

fmc = _^j 1 dx ^^ [1 _ y _^ (615) 

J M<j Ivl c 

with 

Mg D = Ms +^(1 - x )ll ■ (6-16) 
The regularized convergent part of the two-fold ILIs is evaluated by 

r (2)RC _ A f 1 , f ul 2, H+Mf 



-vr 4 f 1 dx [ [djfe?], ^ 

Jo J [ 1J (jfe? + x(l 



x) 9 2 + M z 2 ) 2 



[ ln( 1 + 4 f+M? /VM ) " MK + JWjf) ] (6.17) 
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Consider next the case with a + 7 = 2 and (3 = 1. In this case, the general overlapping loop 
integral contains two logarithmically divergent sub- integrals over k\ (07) and k 2 (Pj)- The overall 
integral is quadratically divergent. In the corresponding ILIs, the sub-integral over the momentum-like 
variable k\ {a + 7 = 2) is logarithmically divergent, whereas the sub-integral over k 2 in ILIs becomes 
quadratically divergent (a + (3 + 7 — 2 = 1) as it characterizes the overall divergences. Similar to the 
first case, it is not difficult to carry out the integration concerning the regularized divergent parts and 
yield the result with the following form 

I ( 2 2)R = jWRDj^RD + jWRCj{1)RD + T (2)RC ^ 

where the regularized quadratically divergent one-fold ILI and the regularized convergent one-fold ILI 
are given by 

I^ RD = vr 2 { Ml - A, [ In ^ - lw + 1 + y 2 (^L) }} (6.19) 

I? RC = f 1 dx I [dk% (6.20) 

2 Jo J 1 J ( k\ + x(l - x)ql + Mf ) V ' 

The regularized convergent two-fold ILI is evaluated by 



Jo J {H + x{l - x)q l Q + Mf 



l 

[ ln( 1 + 4 1+m? /^m ) ~ *S\ + Mf) ] (6.21) 

From the above explicit forms, it is easily seen that when taking the CES M c to be infinite, i.e., 
M c — > 00, we then recover the initial divergent behaviors of the ILIs (or the ones of the corresponding 
general overlapping Feynman integrals). We shall be more interested in the possible harmful diver- 
gences contained in the above results at M c — > 00. One notices that the divergences associated with 
possible non-polynomial external momentum dependent coefficient functions could occur only in the 
double divergent terms fW RD jW RD anc : jW RD l)p RD . This is because only the regularized one-fold 
ILIs Iq RD and j^> RD contain the logarithmic term In M 2 i in which the mass factor Ai 2 is in general 
the function of external momenta. It is remarkable to note that the double divergent terms in the 
two cases contain a common one- fold ILI, i.e., I^ RD . This implies that the two harmful divergences 
at M c — > 00 are actually the same. Therefore, we only need introduce one subtraction term to make 
them becoming harmless. 

Consider now the regularization to the subtraction term for the sub-integral over the loop momen- 
tum k± (see eq.(5.17)) 

" A^' M* -8 + I ld% \ q Jl + M ? f - W*?* («*> 
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which is actually factorized and can easily be carried out. For the case with a + 7 = 2 and j3 = 2, we 
have 

I (2)RS = }(l)RD {mo)I (l)RD {6 23) 

and for the case with a + 7 = 2 and f3 = 1, we yield 

/ (2) RS = jf (l)i?D (mo)/ (l)«D (6>24) 

where the explicit subtraction point m is chosen. 
It then becomes obvious that the differences 

T (2)R A2)RS _ j{l)RC T (1)RD f(l)RC T (1)RD T (2)RC ,„ 
j(2)R _ j(2)RS = j(l)RCj(l)RD + j(l)RC j(l)RD + j(2)RC ^ ^ 

contain only harmless divegencies at M c — > 00. Here Iq is the additional convergent function 
arising from the subtraction 

° - 2 /' [in I + (6.27) 
which is independent of the external momenta and vanishes if choosing = /j% o , 

4 1)i?C = 0, for m 2 = f4 o (6.28) 

The above demonstrations for the two interesting cases have provided an explicit verification on 
the factorization theorem (theorem I) and the subtraction theorem (theorem II) for overlapping di- 
vergences, as well as on the harmless divergence theorem (theorem III) and the trivial convergence 
theorem (theorem IV). Practically speaking, it provides an explicit demonstration on the regulariza- 
tion and renormalization prescriptions for two loop graphs. Obviously, such prescriptions can easily 
be generalized to more closed loops. 

To establish the consistency of the new regularization, we shall further verify the generalized Ward 
identities (or gauge invariant conditions) at two and more loop level. As we have shown in section 
II that only the divergent integrals could spoil the generalized Ward identities. On the other hand, 
the generalized Ward identities only require the regularized tensor and scalar type divergent ILIs to 
satisfy a set of consistency conditions. Therefore the verification is equivalent to check whether the 
consistency conditions between the regularized tensor and scalar type divergent ILIs are still preserved 
at two and more loop level. For this purpose, one only needs to apply the theorems V and VI (i.e., 
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the reduction theorem for overlapping tensor type integrals and the relation theorem for the tensor 
and scalar type ILIs) to the regularized tensor and scalar type divergent ILIs. Obviously, the desired 
result can easily be achieved from those two theorems. 

Last but not least, it must be emphasized that the procedure respects unitarity and causality. 
This is because: (i) the usual Feynman parameter method and the newly formulated UV-divergence 
preserving parameter method do not violate unitarity and causality, this process is independent of 
regularization and can be made in any regularization schemes; (ii) the evaluation of ILIs involves just 
the shifts of integrating loop momenta and the convergent integrals over the loop momenta. This 
procedure is well defined and ensured by the translational invariance and the safety of convergent 
integrals; (iii) the regulators act on the whole Feynman integrals of loop graphs rather than on the 
propagators, which distinguishes to some Pauli-Villars inspired regularization schemes. The Feyn- 
man integrals are well defined in the Euclidean momentum space; (iv) unlike the usual Pauli-Villars 
inspired regularization schemes, the mass factors and numbers of the regulators (or regularization 
quanta) in the new regularization are eventually taken to be infinitely large and decouple from the 
theory. No additional singularities appear via the Feynman propagators; (v) also unlike the BPHZ 
subtraction scheme which is based on expanding around an external momentum, the new regulariza- 
tion and renormalization scheme introduces two intrinsic mass scales, i.e., the CES M c and SES fi s , 
to characterize the UV and IR behaviors of the Feynman amplitudes, so that the structure of the 
corresponding amplitudes involving the external momenta is not changed. This can be seen explicitly 
from the evaluation of the vacuum polarization function presented in section II and appendix A; (vi) 
it has been shown that in a given order any new subtraction terms are finite polynomials in external 
momenta, which is also necessary for maintaining unitarity of the theory. In fact, the subtraction 
process in the new regularization scheme can be made in a similar way as the one in the dimensional 
regularization. 



VII. EVALUATION OF ILIS FOR ARBITRARY LOOP GRAPHS 

It has been seen that evaluating ILIs for loop graphs is a crucial step in the new regularization 
method. We shall present in this section a general description on the evaluation of ILIs for arbitrary 
loop graphs. Though the demonstration may concern some tedious formulae, while it must be very 
useful for a practical computation of more loop diagrams. 

We begin with the scalar type overlapping loop integrals involved in the general n-loop diagrams 
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with n an arbitrary number 



, r, r f 



Ii% J d A k n ^ + n n j ^ {k 2 + M 2 )ai p . _ kj + p. .)2 + M 2 j]aij ( 7 - X ) 

with «j > and ctij > 0. For the loop diagrams in which the internal lines are topologically not 
crossed over each other, the loop integrals are corresponding to the case with j = i + 1 

/i /" 1 1 

+ Mir- n y ^ + - + pu + i? + ai^f^ (7,2) 

Here on > and 0^+1 > as if one of them vanishes, the corresponding sub-integral becomes no 
longer an overlapping one. The most divergence of the above overlapping integrals is corresponding 
to the case with on = 1 (i = 1, • • • n) and au+i = 1 (i = 1, • • • n — 1), which is quadratically divergent. 
This can easily be shown from the power counting rule as the momentum dimension of the overlapping 
integral in this case is two (An — In — (2n — 1) = 2). 

The above overlapping integrals can be regarded as the most general ones for n-loop graphs. Here 
we have also supposed that one can always apply the usual Feynman parameter method to all the 
momentum factors ki , so that all propagators concerning ki and the external momenta pi are combined 
into a single one. Similarly for the propagators concerning (ki — kj) (i ^ j). Therefore, the mass factors 
Ai 2 and Mfj are in general the functions of masses m; and external momenta pi (I = 1, 2, • • • ) 

Mj = M\(m\,p\,m%p\,p x -p 2 ,---), M% = M%(m\,p\,m\,p%p x -p 2 , ■ ■ ■ ) (7.3) 

For simplicity, it is good enough to begin with the demonstration on the evaluation of ILIs for three 
loop graphs. Its generalization to more loop graphs is obvious. The procedure and prescription on the 
evaluation of ILIs for arbitrary loop diagrams will be presented after an explicit demonstration. For 
three loop diagrams, the general overlapping loop integral can be expressed as 

j( 3 ) = f d 4 k 3 I d A k 2 f d 4 k l n7 1 —=r n9 \- ^ —„ — (7.4) 
aiai] J J J (k 2 + Ml) a i {kl + Ml) a * {kl + Ml) a z K ' 

1 1 1 

[(A* -k 2 + p 12 )2 + M\ 2 ] a ^ [{h -k 3 + pi 3 ) 2 + Mj 3 ] a ^ [(k 2 -k 3 + p 23 ) 2 + M 2 23 ] a ™ 

As the first step, we apply the usual Feynman parameter method (see Appendix A) to the denom- 
inators containing the loop momentum k\. Taking the Feynman parameter variables a\ = k\ + M 2 , 
a 2 = (k\ — k 2 + p\ 2 ) 2 + M\ 2 and 03 = (h — k 3 + p\ 3 ) 2 + M\ 3 , the sub-integral over k\ becomes 
quadratic after making a momentum shift 



hn -> hp + (xi ~ x 2 )(k 2 ^ - p 12fl ) + x 2 (k 3/x - p 13fl ) 
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The resulting explicit form is given by 



o 



(1 - Xl) 01 " 1 ^! -X 2 ) ai2 -^X 2 



1~.CH3 — 1 



[ k\ + ( Xl - x 2 )(l - x x + x 2 ){k 2 - pu) 2 + x 2 (l - x 2 )(k 3 - pi 3 ) 2 + M 2 X1X2 ]«i+«ia+aia 

1 1 1 , . 

(7.5) 



( k 2 + M 2 2 )<* [(fca - k 3 + p 23 ) 2 + M 2 3 ]°™ ( fc§ + X 2 ) 



with 



A^ 1Xa = (a* - x 2 )M\ 2 + (1 - xi).M 2 + x 2 M 2 3 (7.6) 

Apply again the usual Feynman parameter method to the remaining denominators containing the 
loop momentum k 2 . Taking the Feynman parameter variables a\ = k 2 + M 2 and a 2 = (k 2 — k 3 + 
P23) 2 + A4 2 3 , and making a momentum shift 

we yield 



4 S L = ^T^M 5" 2 ^ " 23 ! / / d% I d% I dxj ' dx 



r(ai)r(ai 2 )r(ai 3 ) r(a 2 )r(a 23 ) 

(1 - Xi) 01 " 1 ^! - X^^X* 13 ' 1 



'2 



[ k\ + (xi - x 2 )(l - xi + x 2 ){k 2 - Pl2 ) 2 + x 2 (l - x 2 )(A:3 - p 13 ) 2 + M 2 X1X2 ]<*1+<*12+«13 

1 , (l-i/i) 02 " 1 ^ 23 " 1 1 

o yi [ k 2 + yi (l- yi )(k 3 -p l3 ) 2 + Ml }^+^(k 2 + M 2 )^ [ " j 

with 

= (l-yi)X 2 + yi ^ 2 3 (7.8) 

We now adopt as the second step the UV-divergence preserving parameter method. Taking a = 
A;f+(xi-x 2 )(l-xi+x 2 )(A; 2 -pi 2 ) 2 +x 2 (l-x 2 )(A;3-pi3) 2 +A4 2 lX2 and b = k^+y^l-y^ih-p^+M^, 
the integral becomes 

^ = ^^^^^) 1 d% 1 ^ 1 ^ i^r 1 ^ 1 

(1 - x^^-^xi - x 2 ) a ^- l x^ 3 -\l - y!) 02 " 1 ^ 23 " 1 ^ 2 ^ 23 " 1 1 



{ k\ + [ui + (Xi - X 2 )(l - Xi + X 2 )]kl + /i 2 3 + UiM^ + M 2 X1X2 }°i+°i2+«i3+"2+a 2 3 (k 2 + M 2 3 ) a * 
with 

/4 = (xi - x 2 )(l - xi + x 2 )[ 1 - {X1 - X2)(1 - Xl + " 2) - ][ Pl2 - yi (fc 3 - P23 ) ] 2 

Ul + (xi - x 2 )(l - Xl + x 2 ) 

+ [x 2 (l -x 2 ) + yi (l -yi) Ul ](k 3 - Pl3 ) 2 (7.9) 
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where we have made the momentum shift once more 

U 1 (x 1 - x 2 )(l - x ± + x 2 ) 

Ml + (Xl — X 2 )(l - Xl + X 2 ) 

In general, after performing the UV-divergence preserving parameter method, one can safely carry 
out the integration over the momentum k\ as ai + ai 2 -|-ai 3 + a 2 + a 23 > 2 for a non-trivial overlapping 
integral. After integration over k±, the integral gets the following form 

j®. = vr ^i" 1 + q " + " i3 +g + ggg - 2) r d % f d % t l dxi r d x 2 r 1 dyi 



r(ai)r(ai 2 )r(ai 3 )r(a 2 )r(a;23) 



o Jo 



CO „.«2+a23 — 1 1 

1 \ 

Ul [ui + (si - x 2 )(l -xt+ X2 )]( Ql+Ql2+ai3 - 2 ) + ( a2+a23 ) + A^|) Q 3 ^ ' -* 

(1 - xi) 01 - 1 ^! - x 2 ) ai2 - 1 x 2 n3 ~ 1 (l - yi)" 2 - 1 !/" 23 " 1 

{ /c 2 + - Vl ) + Ti^Kfcg - p 23 ) 2 + T U1X [ P12 - Vl (k 3 - p 23 )] 2 + fl 2 ui }«i+«l2+« 13 +«2+«23-2 

with 

rp _ (Xl-X 2 )(l-Xl+X 2 ) (Xi-X 2 )(l-Xl+X 2 ) 

" i:r U\ + (xi — x 2 )(l — Xl + X 2 ) Wl + (xi - x 2 )(l - Xl + x 2 ) 

m _ x 2 (l - x 2 ) - (xi - x 2 ){l -xi + x 2 ) yi (l - yi) 
uixy ~ «i + (xi-x 2 )(l-xi + x 2 ) U j 

^ = ^i + (xi-x 2 )(l-xi + x 2 )- M ' 1 + ui + (xi-x 2 )(l-xi+x 2 ) > '' lX2 (7-13) 
where we have factored out the UV-divergence preserving integration over u\. 

It is of interest to notice that the sub- integral over k 2 in eq.(7.10) becomes solely quadratic. We 
can then repeatedly adopt the UV-divergence preserving parameter method. Taking a = k\ + [yi(l — 
2/1 ) + T lxy ](k 3 - p 23 ) 2 + T uix \p 12 - yi(k 3 - p 23 )] 2 + /4 X and b = k\ + A^§, we have 



r (3) 2 r ( Q l + «12 + "13 + «2 + «23 + «3 ~ 2) f 4 / 4 / f 1 / 

^oi,- = t r / vF? ^F? vr7 s r , s r , , d k 3 d k 2 dx x \ dx 2 / d yi 

13 r(ai)r(ai 2 )r(ai 3 )r(a 2 )r(a 2 3)r(a 3 ) J J Jo Jo Jo 



00 ^C2+Q23~l roo 

dui [ui + (xi - x 2 )(i - xi + x 2 )]( a i+ a i 2+Q i3- 2 )+( Q2+a2 3) y dU2 

(1 - xQ^-^xi - xa)" 12 -^ 13 " 1 ^ - j/i)" 2 " 1 ^ 23 " 1 ^ 3 " 1 

{ k\ + [lt 2 + yi(l - yi) + p ui ] fc 2 + U 2 7W 2 + /i^ + p 2 Ul }<*l+<*i2+ai3+a2+a23+a 3 -2 

with 

?£i = [2/1 (1 - 2/i) + P Ul ] v\ 3 + [p 2 2 - 2yipi2 • p 23 ] 

Pu\ = T UlX yi + T U \xy C^'l^) 
where we have made the momentum shift 

, , , [ 2/1 (1 - 2/l) - TuyxVl + T Ul xy] P23^ + T UlX yip 12fl 

M2 + 2/i(l -yi) + p Ul 
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In general, as a.\ + a.\ 2 + Q13 + a 2 + «23 + «3 — 2 > 2, we can safely perform the integration over 
k 2 and obtain the desired three-fold ILIs 

.(3) _ 4 T(ai + ai2 + ai3 + «2 + «23 + «3 ~ 4) 

a ^ " ^ r(a 1 )r(a 12 )r(a 13 )r(a 2 )r(a 23 )r(a 3 ) 

I dx x r dx 2 (l-x 1 ) a ^ 1 (x 1 -x 2 ) a ^- l x^" 1 [ dytil-y!)"*- 1 ^- 1 
Jo Jo Jo 

OO «2+«23~ 1 

Ul [m + (xi - x 2 )(l - xi + x 2 )]^ 1+ai2+avi - 2)+{ - a2+a2z) 

00 1i 2 3_1 (1) 2 2 

^ 2 [ U2 + yi (l- yi ) + p ui ](«i+«i2+ai3-2)+(a 2 +a 2 3-2)+a 3 I <*i<*ij( M ?> + ^«l« 2 ) ( 7 - 17 ) 

where Ia}aij{M\ + Mu 1U2 ) defines the overall one-fold ILIs 

^nay(-^3 +N« 2 ) = J d k 3r k 2 + M 2 + 2 j(a 1+ai2 + ai3 _2)+(a 2 +a 2 3-2)+a3 ^" 18 ' ) 



with 



/4 - - i/i) + pui] {Ml - p| 3 ) + r uix [pf 2 - 2yip 12 • p 23 ] 



(7.19) 



u 2 + yi(l -yx) + Pu\ 

Similar to the two-fold ILIs of two loop graphs discussed in section V, it is useful to introduce the 
momentum-like integrating variables and measures 



/ * /* d ki 1 



r dui= i 



kf = Ui q z , / dm = / ' „ „ (t = l,2) (7.20) 



Correspondingly, one needs to make the replacement for the relevant quantities 

t UiX — ► t~ 2x , r wixj/ — >■ T 1 ^^ 

2 2 



It is noted that in the infinite limit of kf (or ui) , the quantities T? 2 , Tt 2 , A*| 2 an d p| 2 approach 



to zero. Consequently, we have 



(pui) -»• at /c 2 (ui) -> oo (7.22) 
/"l^i (KUJ^O at A; 2 , fc 2 («i,U2)->ao (7.23) 

This ensures the factorization property for overlapping divergences in three loop graphs. Here the 
sub-integrals over the momentum-like variables k\ and k 2 characterize the UV divergent properties 
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of one-loop and two-loop sub-diagrams respectively, and the integral over the loop momentum 
describes the overall divergent property of three loop diagrams. 

From the above demonstration on the evaluation of ILIs for three loop graphs, we can now straight- 
forwardly generalize it to arbitrary loop graphs. The general prescription on the evaluation of ILIs 
may be summarized as follows: 

(1) . As the first step, one shall repeatedly apply the usual Feynman parameter method to the 
denominators containing the same loop momentum according to the given order say k\ , k%, • • • , so 
that one arrives at the integrals in which n-loop Feynman integral contains correspondingly n's de- 
nominators. By making appropriate momentum shifts, each denominator is governed by the quadratic 
of momentum, say k\ , k^, • • • ■ With this procedure, it is not difficult to see that the momentum in- 
tegral corresponding to the first denominator say k\ becomes purely quadratic. In general, for n-loop 
overlapping integrals, one needs to repeatedly use the Feynman parameter method by n — 1 times. 

(2) . After the step (1), one then adopts the UV-divergence preserving parameter method to the 
first two denominators, i.e., the ones involving k\ and k\. after that, one can safely carry out the 
integration over k\ as the resulting sub- integral in such a way becomes convergent. Note that the 
sub-integral over k\ is quadratic and its integration is easily performed. 

(3) . One first makes momentum shift again for the momentum involved in the second denominator, 
i.e., fe, so that the integral over becomes purely quadratic. Factoring out the coefficient of the 
momentum k%, one arrives at the integral which contains UV-divergence preserving momentum-like 
sub-integral over k\ instead of the one over the loop momentum ki, and the remaining integral has 
the same form as the one resulting from the first step (1) but with the reduced integrals over the loop 
momenta &2, &3> 

(4) . Repeating the steps (2) and (3) to the remaining loop momentum integrals over &2, • • • , 
but keeping the integral over the last loop momentum, we then arrive at the desired ILIs. In general, 
for n-loop overlapping integrals, the resulting ILIs are the n-fold ones with (n-1) sub-integrals over the 
momentum- like variables and one sub- integral over the loop momentum (say k n ). 

According to this procedure and prescription, the resulting ILIs for any loop graphs possess all the 
properties observed in the ILIs of two loop graphs (see section V). In here for the ILIs of more loop 
graphs, one only needs to notice the essential property that all the sub- integrals over the momentum- 
like variables ki, k%^ ■ ■ ■ characterize the divergences of one- loop, two-loop, • • • sub-diagrams, respec- 
tively, and the single sub-integral over the initial loop momentum characterizes the overall divergences 
of the loop graph. One can check that this procedure also ensures that all the UV divergences are 
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characterized by the integrals over the momentum-like variables k\ , fa, • • • , the usual Feynman param- 
eter integrations contain no UV divergences and only IR divergences could be hidden in the Feynman 
parameter integrations. For the regularized ILIs in the new regularization scheme, all the Feynman 
parameter integrations become convergent. It is then obvious that all the theorems deduced in section 
V and the prescriptions for the regularization and renormalization described in section VI hold for 
ILIs of arbitrary loop graphs. 

It is of interest to observe that Feynman integrals of arbitrary loop graphs can always be evaluated 
into the corresponding ILIs which contains only a single one-fold ILI over the initial loop momentum 
which characterizes the overall divergences of the loop graph. Similarly, one can show that the sub- 
tracting the divergences from all the sub- integrals over the momentum-like variables k± , fa, • • • , the 
overall divergence becomes harmless. 



VIII. CONFORMAL SCALING SYMMETRY BREAKING AND 
THE MASS GAP/QUARK CONFINEMENT 

It is particularly noted that the new regularization method may cause the conformal scaling sym- 
metry to be broken down due to the existence of two intrinsic mass scales M c and n s . For instance, 
because of the quadratic 'cut-off' terms of the CES M c , the conformal scaling symmetry can be broken 
down in a class of theories with scalar interactions. In general, as long as taking the SES fj, s to be finite 
(Us 7^ 0), the conformal scaling symmetry will be broken down even in gauge theories. Of particular, 
once a mass gap fi s = fx c is generated via dynamical reasons of strong interactions, it becomes natural 
to set the SES fj, s to be at the mass gap fj, s = fj, c . 

An interesting example is for QCD in which an energy scale around 1 GeV has often be introduced 
to characterize the low energy dynamics of QCD. This is because when the SES fi s runs down to such 
a low energy scale, fi s ~ 1 GeV, the interactions become much stronger. This feature may easily be 
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331 ] . By evaluating the relevant one 



understood from the behavior of asymptotic freedom of QCD 
loop Feynman diagrams, it is not difficult to check that the renormalized gauge coupling constant gn 
from the new regularization method is consistent with the one from the dimensional regularization at 
D — > 4 and M c — > oo 

» w ^ ^ 1 + ^ (&> - f (- f - - + »#) 1 ( "> 

where only the polynomial of the finite terms fj%/M% differ from the one in the dimensional regular- 
ization. In obtaining the above result, the masses of quarks have been ignored. 
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On the other hand, for a finite SES /i s , the conformal scaling symmetry is explicitly broken down 
due to higher dimensional effective interaction terms, i.e., their momentum power counting dimensions 
are larger than four. In general, the UV convergent but IR divergent Feynman diagrams will result in 
such higher dimensional interaction terms. From the dimensional analysis, these higher dimensional 
interactions can be written down in terms of the expansion of the inverse power of the SES Hs 

\{qLq R ){q R qL) , \{G^f ■■■ (8.2) 
K MS 

which implies that the existence of a non-trivial solution for the SES Hs is transmuted to the one of 
nonzero vacuum expectation values of (qq) and (a s G^ u G aiJ,u ). Indeed, it has been turned out that at 
a low energy scale Hs < m c — 1.3 GeV (m c is the charm quark mass), the light quarks and gluons 



were found to have nonzero condensates 



341 ] due to strong interactions 



(qq) ~ (-230MeV) 3 , qo^G^q ~ -(0.9GeV) 
^-(a s G a uv G afM/ ) ~ (238MeV) 4 (8.3) 

47T p 

which implies that both conformal scaling symmetry and chiral symmetry are dynamically broken 
down at such a low energy scale. This strongly indicates that a mass gap with Hs = He ^ m c must 
be generated and the quarks are going to be confined around such a low energy scale. Obviously, the 
conformal symmetry breaking scale or the mass gap Hs = He must be related to the dynamically chiral 
symmetry breaking scale A x which has been known to be around 1 GeV. The order of magnitude for 
the mass gap He may be estimated from the quark and gluon condensates. From simple dimensional 
and large iV c considerations, we have 

(a s G a uu G a ^) 

He ~ - 3 ^ g) ; * 1-lGeV (8.4) 

which holds in the large N c limit with a s N c being held fixed. As expected, the mass gap He an d the 
chiral symmetry breaking scale A x are at the same order of magnitude, namely: 

H s = ^ ~ A x ~ LlGeV. (8.5) 

Correspondingly, the critical temperature for the quark deconfinement is estimated to be at the order 
of magnitude 

T c ~ ^ ~^2L ~ 175MeV. (8.6) 
2vr 2vr K ' 
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Note that to obtain more definitive results, a detailed quantitative calculation is needed. Only the 
perturbative evaluation of QCD is not enough to determine the numerical value for the mass gap. It 
is necessary to develop low energy dynamics of QCD for a nonperturbation evaluation of QCD. One 
may construct some effective perturbation theories to describe the low energy QCD. The crucial issue 
concerns a consistent matching between the perturbative and nonperturbative QCD. In general, the 
IR cut-off in the perturbative QCD is related to the UV cut-off in the nonperturbative QCD. As a 
good example, it has been shown in ref. 



22] that the UV cut-off in the chiral perturbation theory of 
mesons does match to the IR cut-off in perturbative QCD via some appropriate matching conditions. 
Therefore, it is of interest to further investigate the IR and UV correspondence/matching of the 
short and long distance dynamics of QCD. As the new regularization is applicable to both underlying 
perturbative QCD and effective QFTs characterizing nonperturbative QCD at low energy, it may 
provide us a practical tool to understand well the genesis of mass gap and the quark confinement via 
the IR and UV correspondence/matching of the short and long distance dynamics of QCD. 



IX. CONCLUSIONS 



We have explicitly presented a proof for the existence of a symmetry principle preserving and 
infinity free regularization and renormalization method. The main point for the new method is to 
regularize the whole Feynman integrals of loop graphs rather than just the propagators of fields. The 
other important point is to analyze an infinitely large number of regulators (or regularization quanta) 
rather than only a limited few regulators. Specifically, the concepts of irreducible loop integrals (ILIs) 
and UV-divergence preserving parameter method as well as intrinsic mass scales (i.e., the characteristic 
energy scale (CES) M c and the sliding energy scale (SES) fi s ) have been found to play an essential 
role. Consequently, the regularized QFTs become well defined. The method has been developed to 
be practically applicable to both underlying renormalizable QFTs and effective QFTs. Indeed, it 
provides an explicit demonstration that a quantum theory of fields can be applied to describe the 
laws of nature only when the considered energy scale is sufficiently lower than the CES M c . Thus the 
underlying renormalizable QFTs, like Yang-Mills gauge theories, become particularly interesting as 
the CES M c in such theories is in principle allowed to be infinitely large M c — > oo through a suitable 
renormalization for the gauge coupling and fields, which is equivalent to the dimensional regularization 
with D -» 4. 

In general, the scales M c and /i s at a particularly interesting energy scale fj, s = /j, c actually set the 
UV and IR cut-off energy scales respectively in the regularized QFTs, but they distinguish from the 
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naive cutoff momentum scales imposed kinematically to the upper and lower bounds of the internal 
loop momenta, the CES M c and the SES fj, s are more likely to behave as the dynamical mass scales, 
which remains to be understood from a deeper reason. It is intriguing to see that the new regularization 
method provides a consistent way alternative to string theories for making QFTs to be well defined 
finite theories. Though the explicit demonstration for it has been carried out only at one and two loop 
level, its proof to more loops has been shown to be a straightforward generalization. In this sense, the 
divergence problem of QFTs may no longer be a strong reason for extending the underlying theory to 
go beyond the quantum theory of fields /elementary particles. 

First of all, it must be very helpful to further study the possible interesting phenomena caused 
by the conformal scaling symmetry breaking. In particular, the IR and UV correspondence between 
the short distance physics and long distance physics may play an important role to explore the 
nonperturbative effects and to understand well the matching between the short and long distance 
physics. We believe that the conformal scaling symmetry and its breaking mechanism may eventually 
guide us to reveal some longstanding puzzles, such as: genesis of mass gap, unseen quarks, missing 
symmetries, origin of quark and lepton masses and mixing angles, and also small but nonzero 
cosmological constant. Finally, we hope that the new method described in this paper can widen the 
applications of QFTs. 
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APPENDIX A: VACUUM POLARIZATION FUNCTIONS OF GAUGE FIELDS 



In this appendix, we present a detailed evaluation of the gauge field vacuum polarization function 
from various non-vanishing diagrams. In the evaluation, we shall often use the Feynman parameter 
formula 

i _r {a+m f' f-r^ 1 (A1) 



al4 r(a)r(/3)7 [a 1 (l-i) + o 2 x]"+< 3 



and its generalized form 



i r(ai + ••• + «„) f 1 , [ X1 , [ Xn - 2 , 

axi / ax 2 ■ ■ ■ dx n -i 



a^a^---at V{a x ) ■ ■ ■ T(a n ) 



o Jo Jo 



[l-X 1 ) ai - 1 (x 1 -X 2 



\ a '2~ 1 . . . x° n ~^ 



n-1 



(A2) 



[ oi(l - xi) + a 2 (xi - x 2 ) + • • • + a n x n _! 
to simplify the integrals. 

There are in general four Feynman diagrams (see Figures (l)-(4)). Consider first the diagram (1) 
which comes from the Yang-Mills trilinear interaction term. With the gauge boson propagator and 
the Yang-Mills trilinear vertex of the Feynman rules, one has 

TT(l)afe _ 1 2 f f f d k FjilAP, fe) / A q\ 

n uv " "n9 JacdJbdc / TTTTT 77 , , no , ■ w,o | ■ n v A ^J 



2 » JacaJoac J (2vr) 4 ((p + A;) 2 +ie)(A; 2 + z£) 
with 

P^u = IQk^ky + hip^K + Puk^) - 2p fl p v + (hp 2 + 2p ■ k + 2k 2 )g l _ tu 

+ ^[ (k 2 + 2p ■ k - p 2 )k^k u + (k 2 + 3p • k) {pfj,k v + p^) - k% Pu - (k 2 + 2p • k) 2 g^ v ] 

+ — 2 [ (A; 2 - 2p 2 )k^K + (p ■ k){p^k u + p„fc M ) + (p 2 - 2k 2 ) PflPv - (k 2 - p 2 ) 2 g^] 



(p + k) 
X 2 

+ ( p + k ^2 k 2 [ + (P ' k) 2 PnP„ - (p • /c)p 2 (p M fc„ +p u k^) (A4) 



After using the Feynman parameter method with ai = fe and a 2 = (p + fc) , and replacing fc^ by 
(/c M — xp M ), as well as performing some algebra, the gauge field vacuum polarization function Hf^ ab 
can be expressed into the following general form in terms of the ILIs 

nW afe = g 2 d5 ab { [ dx [ I 2 g^ + 5J 2m , + ( (5/2 - 2x{l - x)p 2 g^ - (1 + 5x(l - x)) Py p v ) I } 
Jo 

+ XT (3) dxx[ -hg^v + hp* + ( 2(1 - x)(l + 2x)p 2 g^ + (x 2 - 2)p ii p u ) I 
Jo 

- ( (x + 2(1 - x 2 ) )p 2 g^gl + 4xVp° 9flv - (2x 2 + l)pP (g°p„ + <£p M )) I 0pcr 

■ ( (1 - x) 2 (l + 2x) 2 p A g^ - (1 + x(l - 2x)(2 - x 2 ) )p 2 p ii p v ) J_ 2 ] 

+ A 2 r(4)/2 / x(l - x) [ p 4 g^ + p^p u p p p° - p 2 f (g° Pl/ + fi P(t ) /_ 2pCT ] } (A5) 
J o 
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where C\ and A are defined as 

facdfbcd = Cid a b, A = 1 — £ (A6) 

and the mass term in all the ILIs is given by 

M 2 = -x(l - x)p 2 (A7) 

The diagram (2) is the so-called tadpole diagram with no momentum flowing into the loop and 
it arises from the Yang-Mills quadrilinear interaction terms. This diagram is quadratically divergent 
from power counting and vanishes in dimensional regularization, which has been thought to be the 
crucial point for preserving the gauge invariance in dimension regularization. As a general evaluation 
before applying for any regularization, it is better to treat it in the equal foot without worrying about 
quadratic divergences 



Y[{ 2 ) ab — n 2 f f f 
L1 fiu IJ Jacdjbdc I 



d 4 k 3guv - Kg^v - k^k u /k 2 ) 



(2vr) 4 (k 2 + ie) 



2n x I d 4 k (P + k) 2 [ 2>g^ v - \{g^ y - k^k u /k 2 ) 



- 9 J (2vr)4 {{p - k)2 - i£){k2 - 1£) (A8) 

In order to be able to express this loop integral by using the same ILIs as those appearing in dia- 
gram (1), we have inserted the momentum factor (p + k) 2 /(p + k) 2 . Adopting the same Feynman 
parametrization and performing some algebra, we yield 

ng afe = -g 2 C x b ah { 3 / dx [ h + (1 - x)(l - 2x)p 2 I } 9llu 

Jo 

f 1 1 

+ Ar(3) / dx [ --I 2 9fiV + (1 - x)I 2flv ] 
Jo 1 

+ Ar(3) J dx [ ( ~(1 - x)(l - 2x)p 2 g^ u + (1 - x)x 2 p llPiy J I 
+ ( (1 - x) 2 (l - 2x)p 2 gP l gZ - 2(1 - xfxpf (g°p„ + ) hp* } 

+ XT (3) / dx (I- x) 2 {l - 2x)x 2 P/1 p u p 2 I_2 } (A9) 
Jo 

which shows that without imposing any regularization, the tadpole graphs of Yang-Mills fields are in 
general no vanishing and has actually all divergent structures given by the ILIs. 

We now consider the diagram (3) which arises from the ghost-gauge field interaction term. Using 
the Faddeev-Popov ghost propagator and ghost-gauge field vertex with including a minus sign for the 
closed loop of Grassmann fields, one has 

n( 3 ) a6 — n 2 f f f ^ (ff + k)^k u 

IV " 9 Jacdjbdc J ^ (27r)4 ((p + k)2 + ie)(Jfe2 + . £) 

= -g 2 C 1 5 ab / dx [ I 2 ^ u - x{l - x)p ll p u I } (AlO) 
Jo 
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In obtaining the second equality, we have used the same Feynman parametrization as the previous 
ones. 

Consider finally the diagram (4) due to fermion-gauge field interactions. Using the fermion prop- 
agator and fermion-gauge field vertex with including a minus sign for the closed fermion loop, one 
reads 

d A k (p + k)^k u + (p + k)vkfj, + (m 2 - k 2 - p ■ tyg^ 



= -g 2 ±N f C 2 5ab / dx [ 21 2 ^{m) - h(jn)g yw + 2x(l - x){p 2 g^ u - p^p u )h{m) } (All) 
J o 

where we have used the definition for the group theory factor 

trT a T b = C 2 8 ab (A12) 

Note also that the mass factor in the ILIs is modified by the fermion mass m as indicated in the 
notation If(m) 

M 2 =m 2 - x(l - x)p 2 (A13) 
To simplify the integrals, the following identities are found to be helpful 

j dx x 2 I[x(l — x)] = j dx x(l — x) I[x(l — x)] + - j dx I[x(l — x)] 

JO JO 2 Jq 

/ dx x 3 I[x(l - x)] = ^ / dx x(l - x) I[x{l - x)] + - / dx I[x(l - x)] 
Jo 2 J 2 j 

/ dx x 4 I[x(l - x)} = / dx x 2 (l-x) 2 I[x(l-x)] 
Jo Jo 

+2 / dx x(l - x) I[x(l -x)] + - [ dx I[x(l - x)] (A14) 
Jo 2 Jo 

/ dx x 5 I[x(l — x)] = — j dx x 2 (l — x) 2 7[x(l — x)] 
Jo 2 Jo 

5 f 1 1 f 1 

H — / dx x(l — x) I[x(l — x)] -| — / dx I[x(l — x)] 
2 Jo 2 Jo 



52 



APPENDIX B: SOME REGULARIZED ONE-FOLD ILIS IN CUT-OFF AND DIMEN- 
SIONAL REGULARIZATIONS 



We present here some regularized one-fold ILIs in the cut-off and dimensional regualrizations. The 
relevant one-fold ILIs in the text are 

/d 4 k 1 f d 4 k 1 

(2^4 k 2 -M 2 Io = J (2^0* (k 2 -M 2 ) 2 (B1) 

f d k k^ky f d k k^k v (-an\ 

2 ^ ~ J (2^)4 (k 2 -M 2 ) 2 ° " J Wf (k 2 -M 2 f [ ' 

Consider first the cut-off regularization, the ILIs can be easily evaluated in four dimensional Eu- 
clidean space of momentum by performing a Wick rotation 

p [ K * d 4 k 1 -i . a2 , , A 2 + M 2 . 

fl . /" A2 1 _J_ A 2 + M 2 M 2 

°~ % J (2vr) 4 (-fc2 _ M 2 ) 2 ~ 16vr 2 [ JW 2 + A 2 + M 2 J ( j 

and 

r /" A2 "M* 1 -t ,l l5 1 1/S , >2l A 2 + .M 2 

^ = •7 (2^ (-tf - ^ = 2^ 15? ' 2 A + 4^ ~ ^ ln ' < B5 » 

riJ . f A2 d 4 k -k^k u 1 i r , A 2 + 7W 2 2M 2 M 4 Q 



V (27r) 4 (-A; 2 - M 2 ) 3 4 9flu IGtt 2 M 2 + A 2 + M 2 2(A 2 + M 2 ) 2 

where A 2 is the cut-off momentum. Note that the above integrals are carried out in the Euclidean 
space and the final results are given by rotating back to the Minkowski space, all integrations over 
momentum will be performed in this way except with a specific mention. 

Consider next the dimensional regularization in which the space-time dimension D is taken D = 
4-e 



and 



jR , f d D k -k,K l -» r ,,2^/2-1 r(i + f)r(i-f) 

2 ^ " i (2vr)^ (-fc 2 -M 2 ) 2 " Z? ^ (4vr) B / 2 1 j r(f) 1 j 

* /• d D k -k,k u l 2 ^/2- 2 r(i + f)r(2-f) 

°^ " 7 (2vr) B (-fc 2 -M 2 )3 " D ^ (4tt)^/ 2 ^ j r(f)r(3) 1 J 
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APPENDIX C: MATHEMATICALLY INTERESTING FUNCTIONS AND 
LIMITING NUMBERS INVOLVED IN THE NEW METHOD 

We first present in the following table some numerical results for the functions h w (N), "y w (N) and 
L n (N) when the regulator number N is taken to be typically large. 

Table 1. The quantities h w (N), ^(N) and L n (N) (n = 1,2,3,4,5,6) as functions of the reg- 
ulator number N . 



N 


2 


4 


6 


10 


100 


500 


5000 


10000 


20000 


h w (N) 


1.0407 


1.0614 


1.0676 


1.0715 


1.0667 


1.0586 


1.0487 


1.0462 


1.0439 




1.0198 


0.7645 


0.7059 


0.6640 


0.6037 


0.5927 


0.5860 


0.5848 


0.5839 


Li(N) 


2.079 


1.416 


1.281 


1.187 


1.0560 


1.0325 


1.0182 


1.0158 


1.0138 


L 2 (N) 


6.726 


2.662 


2.048 


1.664 


1.1828 


1.1040 


1.0574 


1.0495 


1.0431 


L 3 (N) 


29.97 


6.37 


4.00 


2.74 


1.4124 


1.2266 


1.1220 


1.1047 


1.0909 


L 4 (N) 


171.74 


18.79 


9.42 


5.27 


1.8072 


1.4218 


1.2191 


1.1868 


1.1613 


HN) 


1209.62 


66.59 


26.24 


11.74 


2.4888 


1.7265 


1.3599 


1.3040 


1.2606 


L 6 (N) 


10141.19 


277.45 


85.08 


30.08 


3.7152 


2.2060 


1.5614 


1.4690 


1.3981 



As the functions behave well. They are expected to approach some finite values at N — > oo, one 
may use a functional fitting method to extract their values at N — > oo with the required precision. 
Here we choose the following fitting function 



Here CT' are the fitting coefficients. In general, one can choose, instead of 1/lniV, other functions 
with vanishing limit at N -> oo. Taking F^N) = h w (N), F 2 (N) = J W (N) and F 2+n (N) = L n (N) 
(n = 1, 2 • • • ), one then has -Fi(oo) = h w , F 2 (oo) = "y w and F 2+n (oo) = L n (n = 1, 2 • • • ). Practically, 
by computing a set of numerical values of the functions h w (N), j w (N) and L n (N) (re = 1, 2 • • • ) for 
some sufficiently large values of N, and solving the corresponding linear equations, one is then able to 




n=l 



(CI) 
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obtain F a (oo) with the required precision. Here we only present the solutions at the precision 10 



h w = -Fi(oo) ~ 1.000 (C2) 

7™ = F 2 (oo) ~ 5.772 x 10 _1 (C3) 

Li = F 3 (oo) ~ 1.000 (C4) 

L 2 = F 4 (oo) ~ 1.000 (C5) 

L 3 = F 5 (oo) ~ 1.000 (C6) 

L 4 = F 6 (oo) ~ 1.000 (C7) 

L 5 = F 7 (oo) ~ 1.000 (C8) 

L 6 = F 8 (oo) ~ 1.000 (C9) 

They are resulted from a set of values for the functions h w (N), "f w (N) and L n {N) (n = 1,2- ••) 



obtained for N from N = 300 to N = 10000. For a consistent check, such solutions are found to 
provide, with the same precision, a prediction for the numerical values of the functions h w (N), j w (N) 
and L n (N) (n = 1, 2 • • • ) at N = 20000 (see table). 

The above numerical analysis provides a reasonable check on the three conjectures made in section 

IV. 

APPENDIX D: SOME USEFUL REGULARIZED ONE-FOLD ILIS 
IN THE NEW METHOD 

Using the definitions and prescriptions of the new regularization, we present the following useful 
regularized ILIs 
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with 



where 



/ 



1 M 2 u 2 



d 4 k] 



d A k]i 



d 4 k]i 



1 



M 



ln ^-7- + yo(^o) } 



{k 2 + Mf + M 2 f 
1 



A* 



M 2 C 



1 2 M 



1 



(k 2 + Mf + X 2 ) c 



7T 



2 r(a-2) 1 



l-J/-2(a-2)(xrt) } 



r(a) (a* 2 ) 



M?' 



a > 3 



k fjjkp 



M 2 



^ 2 {M 2 -^[ln^- lw + l + y 2 (^) }} 



(k 2 + Mf + X 2 ) 2 2 



M 2 ' 



kfj/k-i/ 

" 1 {k 2 + Mf + X 2 ) 3 ~ 4 



1 r 2 r , M c /Mm 

7 ff^7r { In - 7«- + ?/o(-^2 ) } 



d 4 £^ 



k^xk^ 



11 fji Z 
( fc 2 + M 2 + ^2)4 = 4 Vg^T 2 { 1 " ^(j^) > 



k jjbki, 



1 „ 2 r(a-2) 1 



(A; 2 + Mf + .M 2 ) 



l-|/-2(a-2)(vW) } 



r(a + l) (m 2 ) 



,4 , i knk v knk a 
dk\i 



1 - 2ru2 2r i_ M c , -i , .. / A* 



(fc 2 + Mf + A^ 2 )3 " 8^^ {Mc " M In ^ - 7. + 1 + 2/2(^) ]} 



M 2 ' 

2 

M?- 



^4^j kiik u kpk a 



M 2 

(k 2 + Mf + = 24 WF 2 { ln ^# - 7» + Ito(^) } 



A* 



(k 2 + M 2 + M 2 f 2A" { » vpa} An 2 " L ^ y " zv M 2 



^ 2 {l-y-2(^)l 



d 4 £^ 



/c /l^/ kp kfj 



—6. 



(k 2 + Mf + 7W 2 ) a + 2 4 



7T 



,r(a-2) 1 
r(a + 2) (/X 2 )"" 2 



{ l-y-2(a-2)(^2) } 



(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 

(D10) 

(Dll) 

(D12) 
(D13) 



N 

[d**],= lim £cf 

1=0 



[ d A k= lim Y(-l) l rAT Nl ni [ 



d 4 k 



(D14) 



Mf = /x 2 + ZM| , Ml = M 2 h w (N) In iV 
fx 2 = /x 2 + .M 2 

.. f w ( x/{lh w {N)\uN) \ 
l/-2( Q -2) (x) = - hm c t { 1 + x/{lhw(N)lnN) ) 



i=i 



a-2 



1 W = 1E = 0.5772 • 



h w (N -»■ oo) = 1 



(D15) 
(D16) 
(D17) 

(D18) 
(D19) 



5G 



The explicit forms for y2(%), Vo( x ) an d y-2(x) are given in eqs.(4.10-4.13). 
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Figures 
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(2) (4) 



Vacuum polarization diagrams of gauge field: (1) from trilincar gauge boson interactions ; 
(2) from quartic gauge boson interactions; (3) from gauge-ghost boson interactions; 
(4) from fermion-gaugc boson interactions. 
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